Cell-free nucleic acids for the analysis of the human microbiome and components thereof

ABSTRACT

Methods, devices, compositions and kits are provided for analysis of the microbiome or individual components thereof in an individual. The methods find use in a determination of infection, in analysis of the microbiome structure, in determining the immunocompetence of an individual, and the like. In some embodiments of the invention, the individual is treated with an therapeutic regimen, e.g. drugs, diet, radiation therapy, and the like.

CROSS REFERENCE

This application claims benefit and is a Continuation of applicationSer. No. 16/554,286, filed Aug. 28, 2019, which claims benefit ofapplication Ser. No. 14/536,300 filed Nov. 7, 2014, which claims benefitof U.S. Provisional Patent Application No. 61/901,114 filed Nov. 7,2013, and U.S. Provisional Patent Application No. 61/901,857, filed Nov.8, 2013, which applications are incorporated herein by reference intheir entirety.

GOVERNMENT RIGHTS

This invention was made with Government support under grant RC4AI092673awarded by the National Institutes of Health. The Government has certainrights in this invention.

BACKGROUND OF THE INVENTION

The human microbiome is now recognized as an important component ofhuman health. Community level analyses have shed light on factors thatshape the structure of the bacterial and viral components of themicrobiome, such as age, diet, geographical location, antibiotictreatment and disease. For example, an individual microbiome can bealtered by infection with a pathogenic organism, such that there is anincreased prevalence of that organism systemically, or in an undesirabletissue. The microbiome can also be altered by changes in theimmunocompetence of the individual.

For a variety of purposes it would be desirable to have a method forrapid identification of specific microbiome components, e.g. thepresence and prevalence of commensal, mutualistic, parasitic,opportunistic and pathogenic organisms in an individual microbiome; aswell as an analysis of the overall microbiome structure. The presentinvention provides sensitive, rapid, non-invasive methods of monitoringorgan transplant patients for immunocompetence that are sensitive, rapidand inexpensive.

SUMMARY OF THE INVENTION

The invention provides methods, devices, compositions and kits foranalysis of the microbiome or individual components thereof in anindividual. The methods find use in a determination of infection, inanalysis of the microbiome structure, in determining theimmunocompetence of an individual, and the like. In some embodiments,the invention provides methods of determining the presence andprevalence of microorganisms in an individual, comprising the steps of:(i) providing a sample of cell-free nucleic acids, i.e. DNA and/or RNAfrom an individual; (ii) performing high-throughput sequencing, forexample from about 10⁵, and up to about 10⁹ or more reads per sample;(iii) performing bioinformatics analysis to subtract host sequences,i.e. human, cat, dog, etc. from the analysis; and (iv) determining thepresence and prevalence of microbial sequences, for example by acomparison of the coverage of sequences mapping to a microbial referencesequence to coverage of the host reference sequence.

The subtraction of host sequences may include the step of identifying areference host sequence, and masking microbial sequences or microbialmimicking sequences present in the reference host genome. Similarly,determining the presence of a microbial sequence by comparison to amicrobial reference sequence may include the step of identifying areference microbial sequence, and masking host sequences or hostmimicking sequences present in the reference microbial genome.

A feature of the invention is the unbiased analysis of cell-free nucleicacids from an individual. The methods of the invention generally includean unbiased amplification step, for example by performing PCR withuniversal primers, or by ligation of adapters to the nucleic acid andamplifying with primers specific for the adaptors. The methods of theinvention are typically performed in the absence of sequence specificamplification of microbial sequences. A benefit to this approach is thatanalysis then includes all available microbiome sequences, however itrequires bioinformatics analysis to identify sequences of interest in acomplex dataset predominated by host sequences.

A further benefit of the methods of the invention is the ability toprovide a rapid assessment of an individual microbiome, for exampleanalysis may be completed in less than about 3 days, less than about 2days, less than one day, e.g. less than about 24 hours, less than about20 hours, less than about 18 hours, less than about 14 hours, less thanabout 12 hours.

The analysis of the presence and prevalence of microbial sequences canbe used to provide a determination on infection, of response to therapy,including anti-microbial treatment such as treatment with antibiotics,anti-viral agents, immunization, passive immunotherapy, and the like;diet; immunosuppression, and the like; of response in clinical trials,etc. The information obtained from the analysis may be used to diagnosea condition, to monitor treatment, to select or modify therapeuticregimens, and to optimize therapy. With this approach, therapeuticand/or diagnostic regimens can be individualized and tailored accordingto the specificity data obtained at different times over the course oftreatment, thereby providing a regimen that is individually appropriate.In addition, patient samples can be obtained at any point during thetreatment process; following exposure to a pathogen; course ofinfection, etc. for analysis. The analysis of presence and prevalence ofmicrobial sequences can be provided as a report. The report may beprovided to the individual, to a health care professional, etc.

In some embodiments, the cell-free nucleic acid is obtained from abiological sample selected from the group consisting of blood, serum,cerebrospinal fluid, synovial fluid, urine, and stool. The nucleic acidis extracted from the cell free portion of the sample, e.g. serum orplasma portion of blood may be used. In some embodiments, the nucleicacid is selected from the group consisting of double-stranded DNA,single-stranded DNA, single-stranded DNA hairpins, DNA/RNA hybrids, RNAand RNA hairpins. In some embodiments, the nucleic acid is selected fromthe group consisting of double-stranded DNA, single-stranded DNA andcDNA. In some embodiments, the nucleic acid is mRNA. In someembodiments, the nucleic acid is circulating cell-free DNA.

In some embodiments, the methods comprise quantitating the one or morenucleic acids to provide a determination of prevalence of amicroorganism in the sample. In some embodiments, the amount of the oneor more nucleic acids above a predetermined threshold value isindicative of infection or altered prevalence. In some embodiments,there are different predetermined threshold values for differentmicrobial organisms. In some embodiments, temporal differences in theamount of the one or more nucleic acids are indicative of changes ininfection, altered prevalence, response to therapy, etc.

In some embodiments, the invention provides computer readable mediumscomprising: a set of instructions recorded thereon to cause a computerto perform the steps of: (i) receiving high throughput sequencing datafrom one or more nucleic acids detected in a sample of cell-free nucleicacids from a subject; (ii) performing bioinformatics analysis tosubtract host sequences, i.e. human, cat, dog, etc. from the analysis;and (iii) determining the presence and prevalence of microbialsequences, for example by a comparison of the coverage of sequencesmapping to a microbial reference sequence to coverage of the hostreference sequence.

In some embodiments, the invention provides reagents and kits thereoffor practicing one or more of the methods described herein.

In some embodiments, compositions and method are provided for theassessment of immunocompetence of an individual, particularly anindividual human, by analysis of the microbiome, for example by analysisof the virome. In some embodiments of the invention, the individual istreated with an immunosuppressive regimen, e.g. drugs, radiationtherapy, and the like. In some embodiments the individual is a graftrecipient treated with an immunosuppressive regimen. In some embodimentsthe individual has an autoimmune disease treated with animmunosuppressive regimen. In other embodiments an individual isassessed for immunocompetence in the absence of an immunosuppressiveregimen.

In some embodiments a measurement from an individual is taken at two ortime points, where a change in virus burden is indicative of a change inimmunocompetence. The individual may be treated in accordance with theassessment of immunocompetence, e.g. where an indication of undesirableincreased immunocompetence in a transplant patient is treated withincreased levels of immunosuppressive agents; or where an undesirabledecrease in immunocompetence is treated with therapeutic agents, e.g.anti-viral agents, etc.

Nucleic acid analysis is used to identify and quantify nonhumancell-free nucleic acids in a sample collected from a patient. Thecomposition of the components of the microbiome is performed asdescribed above. The structure of the viral component of the microbiome(the virome) allows a prediction of immunocompetence. In someembodiments, the methods further comprise establishing a virome profileprior to an immunosuppressive regimen, at the initiation of animmunosuppressive regimen, or during the course of an immunosuppressiveregimen, which is used as a reference to changes in the individualvirome. In some embodiments the circulating cell-free DNA isannellovirus DNA.

In particular, the load of viruses of the anelloviridae family is apredictor of immune strength, which is correlated with the probabilityof organ transplant rejection. While other viruses may also bepredictive, it is common for patients to be treated with antivirals thataffect the load of such viruses.

In some embodiments, the invention provides methods of diagnosing orpredicting transplant status or outcome comprising the steps of: (i)providing a sample from a subject who has received a transplant from adonor; (ii) determining the presence or absence of one or more viromenucleic acids; and (iii) diagnosing or predicting transplant status oroutcome based on the virome load. In some embodiments, the transplantstatus or outcome comprises rejection, tolerance, non-rejection basedallograft injury, transplant function, transplant survival, chronictransplant injury, or titer pharmacological immunosuppression. In someembodiments, the amount of the one or more nucleic acids above apredetermined threshold value is indicative of viral load andimmunocompetence. In some embodiments, the threshold is a normativevalue for clinically stable post-transplantation patients with noevidence of transplant rejection or other pathologies. In someembodiments, there are different predetermined threshold values fordifferent transplant outcomes or status. In some embodiments, temporaldifferences in the amount of the one or more nucleic acids areindicative of immunocompetence.

In any of the embodiments described herein, the transplant graft maybeany solid organ, bone marrow or skin transplant. In some embodiments,the transplant is selected from the group consisting of kidneytransplant, heart transplant, liver transplant, pancreas transplant,lung transplant, intestine transplant and skin transplant.

In some embodiments, the invention provides reagents and kits thereoffor practicing one or more of the methods described herein.

All publications and patent applications mentioned in this specificationare herein incorporated by reference to the same extent as if eachindividual publication or patent application was specifically andindividually indicated to be incorporated by reference.

BRIEF DESCRIPTION OF THE DRAWINGS

The novel features of the invention are set forth with particularity inthe appended claims. A better understanding of the features andadvantages of the present invention will be obtained by reference to thefollowing detailed description that sets forth illustrative embodiments,in which the principles of the invention are utilized, and theaccompanying drawings of which:

FIG. 1. Study design, read statistics and phylogenetic distribution. A.Immunosuppression reduces the risk of rejection in transplantation butincreases the risk of infection. B. Design of study. 656 plasma sampleswere collected, cell-free DNA was purified and sequenced to an averagedepth of 1.2 Gbp per sample. C. Number of samples collected as functionof time for the different patient groups part of the study. D. Treatmentprotocol for patients in the study cohort, all patients are treated withmaintenance immunosuppression (tacrolimus-based (TAC) for adult heartand lung transplant recipients and cyclosporine (CYC) for pediatricpatients). CMV positive (donor or recipient, CMV+) transplant cases aretreated with anti-CMV prophylaxis, valganciclovir (VAL). Mean level oftacrolimus measured in blood of transplant recipients treated with aTAC-based protocol (dashed line actual, solid line window averagefilter). E. Fraction of reads that remain after filtering of lowerquality and duplicate reads (mean 86%, left) and after removal of humanand low complexity reads (mean 2%, right). F. Relative genomic abundanceat different levels of taxonomic classification after removal of humanreads (average over all samples from all organ transplant recipients(n=656)).

FIG. 2. Relative viral genomic abundance as a function of drug dose andcomparison to healthy reference. A. Mean virome composition for patientstreated with the immunosuppressant tacrolimus (47 patients, 380 samples)as function of antiviral drug dose (valganciclovir) and concentrationtacrolimus measured in blood. To account for the delayed effect of thevirome composition on drug dose, the data on drug doses were windowaverage filtered (window size 45 days, see FIG. 1C). Herpesvirales andcaudovirales dominate the virome when patients receive low doses ofimmunosuppressants and antiviral drugs. Conversely, anelloviridaedominate the virome when patients receive high doses of these drugs. B.Comparison of virome composition corresponding to healthy references(n=9), post-transplant day one samples with low drug exposure, (n=13),and samples corresponding to high drug exposure (tacrolimus≥9 ng/ml,valganciclovir≥600 mg, n=68). The virome structure for day one samples(1) and the virome structure measured for a set of healthy individuals(H) are distinct from the anellovirus-dominated distribution measuredfor samples corresponding to high drug doses (D). The piecharts show themean fractions, p-values in boxplot based on the Mann-Whitney test. C.Bray-Curtis beta diversity for all samples, among patients with the sametransplant type (heart or lung), within subjects, for patients treatedwith a similar drug dosage (tacrolimus level±0.5 ng/ml,valganciclovir±50 mg), and for samples collected from the same subjectswithin a one-month timespan.

FIG. 3. Temporal dynamics of the microbiome composition post-transplant.A. Relative abundance of dsDNA and ssDNA viruses for different timeperiods (average for all samples). The relative abundance of ssDNAviruses increases rapidly after the onset of the post-transplant drugtherapy. After 6 months, the opposite trend is observed. B. Viral genomeabundance at the family and order level of taxonomic classification fordifferent time periods. The fraction of anelloviridae expands rapidly inthe first several months post-transplant. The fraction of herpesvirales,caudovirales and adenoviridae decreases in that same time period. After6 months, the opposite trends are observed. C. Time-variation in therelative abundance of bacterial phyla. Compared to the viral abundance,the representation of different bacterial phyla is relatively unchangedover the observed post-transplant period. D. Shannon entropy as ameasure of the within-sample alpha-diversity for bacterial and viralgenera as function of time (data grouped per one month time period).

FIG. 4. Virome composition and total viral burden in the absence andpresence of antiviral prophylaxis. A. Absolute viral load as a functionof time, measured as viral genome copies per human genome copiesdetected by sequencing. Box plots are shown for different time periodswith centers of the time periods marked on the x-axis. For all patientclasses, the total viral load increases in the first weekspost-transplant (black line is sigmoid fit, change in load 7.4±3). B.Viral load and composition for CMV+ cases that are treated with bothimmunosuppressants and antiviral drugs (78 patients, 543 samples). C.Viral load and composition for CMV −/− cases, only treated withimmunosuppressants (12 patients, 75 samples).

FIG. 5. Lower anellovirus burden in patients that suffer from graftrejection. A. Time dependence of the anellovirus load in the subgroup ofpatients that suffer from a severe rejection episode (biopsy grade2R/3A, red data, 20 patients, 177 time points) and in the subgroup ofpatients that do not suffer from a severe-rejection-free post-transplantcourse (blue data, 40 patients, 285 time points). Box plots are shownfor different time periods with centers of the time periods marked onthe x-axis. Solid lines are cubic splines (smoothing parameter 0.75).The inset shows a cartoon of the expected opposite association of theincidence of rejection and infection with immunocompetence. B.Anellovirus load relative to the average load measured for all samplesat the same time point. The time-normalized load for non-rejectingpatients (N=208) is compared to the load measured for patients sufferingfrom a mild rejection event (biopsy grade 1R, N=102) and patientssuffering from a severe rejection episode (biopsy grade 2R/3A, N=22).The p-values reflect the probability that the median viral load ishigher for the subgroups at greater risk of rejection. The p-values arecalculated by random sampling of the population with a greater amount ofmeasurement points. N-fold random sampling,p=sum(median(A_(rej))>median(A_(non-rej)))/N), where N=10⁴ and A_(rej)and A_(non-rej) are the relative viral loads for the populations atgreater and lesser risk of rejection and non-rejecting respectively. CTest of the performance of the relative anellovirus load in classifyingpatients as non-rejecting vs. severely rejecting, receiver-operatingcharacteristic curve, area under the curve=0.72.

FIG. 6. Genome Sizes and Hit Statistics, qPCR Assay and Influence ofRead Length on Measured Relative Abundance of Species at DifferentLevels of Taxonomic Classification (A) Distribution of genome sizes inthe reference database with, 1401 viral genomes, 32 fungal genomes and1980 bacterial genomes. (B) Distribution of unique blast hits permillion unique molecules sequenced (average number of hits specified inthe x axis label). C. Distribution of genome equivalents (infectiousagents/diploid human) for viruses, bacteria and fungi (average number ofgenome equivalents specified in the x axis label). (D) Comparison ofsequencing hits found per million of total reads sequenced to the numberof viral copies detected using qPCR. For the qPCR assays, DNA waspurified from 1 ml of plasma and eluted in a 100 μl volume. (E)Measurements of CMV and parvovirus burdens for selected cases. Thehighest loads of CMV virus (genome equivalents, viral/human diploid,G.E.) measured for all samples corresponded to two cases of clinicallydiagnosed disseminated CMV infection (a and b, shaded area denotestime-window of clinical diagnosis, * denotes time of death). (c) shows atime trace of a pediatric patient that suffered from CMV viremia.Parvovirus was detected in one pediatric heart transplant patientimmediately post-transplant (d).* (F) Influence of read length onmeasured relative abundance of species at different levels of taxonomicclassification (n=52). Spearman sample-to-sample correlation, r, and pvalue, p, (two-sample Mann-Whitney U test) for the abundance of the mostabundant node extracted from the 50 and 100 bp data sets: r=0.80, p=0.8(a), r=0.86, p=0.4 (b), r=0.92, p=0.6 (c), r=0.84, p=0.5 (d), r=0.7,p=0.28 (e), r=0.99, p=1 (f).

FIG. 7. Average Drug Doses and Measured Levels for Adult Heart and LungTransplant Patients after Transplantation and Influence of Drug Dosageon Virome composition. (A-C) Average dose of valganciclovir andprednisone (A and C) administered and measured level of tacrolimus inblood (B) for the adult heart and lung transplant patients part of thisstudy. (D) Compared to the viral component, the composition of thebacterial component of the microbiome is relatively insensitive toantivirals and immunosuppressants. (E) Virome composition as function ofdose of anti-CMV drug (valganciclovir) and immunosuppressant(prednisone).

FIG. 8. Temporal Dynamics of the Bacterial Component of the MicrobiomePosttransplant, (A) Relative abundance of bacterial phyla as function oftime. (B) Relative abundance of bacterial genera as function of time.

FIG. 9. Virome Composition and Total Viral Burden for Different PatientClasses, (A and B) Viral load and composition for CMV positive adultheart (A), adult lung (B), and pediatric heart (C) transplantrecipients, treated with both immunosuppressants and antivirals.

FIG. 10. CMV-infection induced allograft injury. A. Correlation betweenclinical report of CMV (Human herpes virus 5, HHV-5) infection fromspecific body fluids (BAL and serum) with donor organ cfdDNA signalmatched to clinical test date (P-values; Mann-Whitney U test). B.P-values for the correlation between clinical diagnosis of infection andcell-free DNA level (dashed line indicates the bonferonni-correctedsignificance threshold) for infections with greater than one clinicalpositive test result. C. An ROC curve that tests the performance of CMVderived cell-free DNA level in CMV-positive and CMV-negative patients(AUC=0.91).

FIG. 11. Monitoring the infectome. A. Clinical testing frequencycompared to the incidence of viral infections detected in sequencing. B.Time-series data for patients who tested positive (red arrows) forspecific infections relative to those who were un-tested. (1) Adenovirussignal in L78 with clinical positives highlighted relative to untestedpatient (L34). (2) Polyomavirus signal in L69 with one positive testrelative to sustained signal in un-tested patient (L57). (3) Threeherpesvirus infections (HHV-4, 5, and 8) in L58 with both positive (red)and negative (black) tests for CMV (HHV-5) highlighted. (4)Microsporidia signal in I6, with four positive tests shown, relative tothe signal observed in L78, which had symptoms of microsporidiosis, butwas un-tested. Data are logged genome equivalents relative to humanwhere zero values were replaced with the detection limit of the assay(the number of genome equivalents consistent with a single sequence readassigned to the target genome).

DETAILED DESCRIPTION OF THE INVENTION

Reference will now be made in detail to particularly preferredembodiments of the invention. Examples of the preferred embodiments areillustrated in the following Examples section.

Unless defined otherwise, all technical and scientific terms used hereinhave the same meaning as is commonly understood by one of skill in theart to which this invention belongs. All patents and publicationsreferred to herein are incorporated by reference in their entirety.

Where a range of values is provided, it is understood that eachintervening value, to the tenth of the unit of the lower limit unlessthe context clearly dictates otherwise, between the upper and lowerlimit of that range and any other stated or intervening value in thatstated range, is encompassed within the invention. The upper and lowerlimits of these smaller ranges may independently be included in thesmaller ranges and are also encompassed within the invention, subject toany specifically excluded limit in the stated range. Where the statedrange includes one or both of the limits, ranges excluding either orboth of those included limits are also included in the invention.

Certain ranges are presented herein with numerical values being precededby the term “about.” The term “about” is used herein to provide literalsupport for the exact number that it precedes, as well as a number thatis near to or approximately the number that the term precedes. Indetermining whether a number is near to or approximately a specificallyrecited number, the near or approximating unrecited number may be anumber which, in the context in which it is presented, provides thesubstantial equivalent of the specifically recited number.

The practice of the present invention employs, unless otherwiseindicated, conventional techniques of immunology, biochemistry,chemistry, molecular biology, microbiology, cell biology, genomics andrecombinant DNA, which are within the skill of the art. See Sambrook,Fritsch and Maniatis, MOLECULAR CLONING: A LABORATORY MANUAL, 2ndedition (1989); CURRENT PROTOCOLS IN MOLECULAR BIOLOGY (F. M. Ausubel,et al. eds., (1987)); the series METHODS IN ENZYMOLOGY (Academic Press,Inc.): PCR 2: A PRACTICAL APPROACH (M. J. MacPherson, B. D. Hames and G.R. Taylor eds. (1995)), Harlow and Lane, eds. (1988) ANTIBODIES, ALABORATORY MANUAL, and ANIMAL CELL CULTURE (R. I. Freshney, ed. (1987)).

The invention provides methods, devices, compositions and kits foranalysis of the microbiome or individual components thereof in anindividual. The methods find use in a determination of infection, inanalysis of the microbiome structure, in determining theimmunocompetence of an individual, and the like. In some embodiments,the invention provides methods of determining whether a patient orsubject is displaying immunocompetence. The term “individual”, “patient”or “subject” as used herein includes humans as well as other mammals.

Definitions.

As used herein the term “diagnose” or “diagnosis” of a status or outcomeincludes predicting or diagnosing the status or outcome, determiningpredisposition to a status or outcome, monitoring treatment of patient,diagnosing a therapeutic response of a patient, and prognosis of statusor outcome, progression, and response to particular treatment.

Microbiota. As used herein, the term microbiota refers to the set ofmicroorganisms present within an individual, usually an individualmammal and more usually a human individual. The microbiota may includepathogenic species; species that constitute the normal flora of onetissue, e.g. skin, oral cavity, etc., but are undesirable in othertissues, e.g. blood, lungs, etc.; commensal organisms found in theabsence of disease; etc. A subset of the microbiome is the virome, whichcomprises the viral components of the microbiome.

The term “microbiome component” as used herein refers to an individualspecies, or in the instance of some gut commensal organisms, a genus ofmicroorganisms. The component may be a viral component, a bacterialcomponent, a fungal component, etc.

In a healthy animal, while the internal tissues, e.g. brain, muscle,etc., are normally presumed to be relatively free of bacterial species,the surface tissues, i.e., skin and mucous membranes, are constantly incontact with environmental organisms and become readily colonized byvarious microbial species. The mixture of organisms known or presumed tobe found in humans at any anatomical site is referred to as the“indigenous microbiota”, including various components of the indigenousmicrobiota. In addition to the indigenous microorganisms are varioustransient components, such as pathogenic or opportunistic infections.Reference sequences of organisms described below are publicly availableand known in the for, example at the Genbank database.

The intestinal microbiota of humans is dominated by species found withintwo bacterial phyla: members of the Bacteroidetes and Firmicutes makeup >90% of the bacterial population. Actinobacteria (e.g., members ofthe Bifidobacterium genus) and Proteobacteria among several other phylaare less prominently represented. Common species of interest includeprominent or less abundant members of this community, and may comprise,without limitation, Bacteroides thetaiotaomicron; Bacteroides caccae;Bacteroides fragilis; Bacteroides melaninogenicus; Bacteroides oralis;Bacteroides uniformis; Lactobacillus; Clostridium perfringens;Clostridium septicum; Clostridium tetani; Bifidobacterium bifidum;Staphylococcus aureus; Enterococcus faecalis; Escherichia coli;Salmonella enteritidis; Klebsiella sp.; Enterobacter sp.; Proteusmirabilis; Pseudomonas aeruginosa; Peptostreptococcus sp.; Peptococcussp., Faecalibacterium sp.; Roseburia sp.; Ruminococcus sp.; Dorea sp.;Alistipes sp.; etc.

In the skin microbiome most bacteria fall into four different phyla:Actinobacteria, Firmicutes, Bacteroidetes and Proteobacteria.Microorganisms that are generally regarded as skin colonizers includecoryneforms of the phylum Actinobacteria (the genera Corynebacterium,Propionibacterium, such as Propionibacterium acnes; and Brevibacterium),the genus Micrococcus and Staphylococcus spp. The most commonly isolatedfungal species are Malassezia spp., which are especially prevalent insebaceous areas. The Demodex mites (such as Demodex folliculorum andDemodex brevis) may also be present. Other types of fungi that arethought to grow on the skin, include Debaryomyces and Cryptococcus spp.As non-commensals, burn wounds commonly become infected with S.pyogenes, Enterococcus spp. or Pseudomonas aeruginosa, and can alsobecome infected with fungi and/or viruses. S. epidermidis is a verycommon skin commensal, but it is also the most frequent cause ofhospital-acquired infection on in-dwelling medical devices such ascatheters or heart valves. For a review, see Nat Rev Microbiol. (2011)April; 9(4):244-53.

Pathogenic species may be bacteria, virus, protozoan parasites, fungalspecies, etc. Bacteria include Brucella sp., Treponema sp.,Mycobacterium sp., Listeria sp., Legionella sp., Helicobacter sp,Streptococcus sp, Neisseria sp, Clostridium sp, Staphylococcus sp. orBacillus sp.; including without limitation Treponema pallidum,Mycobacterium tuberculosis, Mycobacterium leprae, Listeriamonocytogenes, Legionella pneumophila, Helicobacter pylori,Streptococcus pneumoniae, Neisseria meningitis, Clostridium novyi,Clostridium botulinum, Staphylococcus aureus, Bacillus anthracis, etc.

Parasite pathogens include Trichomonas, Toxoplasma, Giardia,Cryptosporidium, Plasmodium, Leishmania, Trypanosoma, Entamoeba,Schistosoma, Filariae, Ascaria, Fasciola; including without limitationTrichomonas vaginalis, Toxoplasma gondii, Giardia intestinalis,Cryptosporidium parva, Plasmodium falciparum, Trypanosoma cruzi,Entamoeba histolytica, Giardia lamblia, Fasciola hepatica, etc.

Viruses that infect humans include, for example, Adeno-associated virus;Aichi virus; Australian bat lyssavirus; BK polyomavirus; Banna virus;Barmah forest virus; Bunyamwera virus; Bunyavirus La Crosse; Bunyavirussnowshoe hare; Cercopithecine herpesvirus; Chandipura virus; Chikungunyavirus; Cosavirus A; Cowpox virus; Coxsackievirus; Crimean-Congohemorrhagic fever virus; Dengue virus; Dhori virus; Dugbe virus;Duvenhage virus; Eastern equine encephalitis virus; Ebolavirus;Echovirus; Encephalomyocarditis virus; Epstein-Barr virus; European batlyssavirus; GB virus C/Hepatitis G virus; Hantaan virus; Hendra virus;Hepatitis A virus; Hepatitis B virus; Hepatitis C virus; Hepatitis Evirus; Hepatitis delta virus; Horsepox virus; Human adenovirus; Humanastrovirus; Human coronavirus; Human cytomegalovirus; Human enterovirus68, 70; Human herpesvirus 1; Human herpesvirus 2; Human herpesvirus 6;Human herpesvirus 7; Human herpesvirus 8; Human immunodeficiency virus;Human papillomavirus 1; Human papillomavirus 2; Human papillomavirus16,18; Human parainfluenza; Human parvovirus B19; Human respiratorysyncytial virus; Human rhinovirus; Human SARS coronavirus; Humanspumaretrovirus; Human T-lymphotropic virus; Human torovirus; InfluenzaA virus; Influenza B virus; Influenza C virus; Isfahan virus; JCpolyomavirus; Japanese encephalitis virus; Junin arenavirus; KIPolyomavirus; Kunjin virus; Lagos bat virus; Lake Victoria marburgvirus;Langat virus; Lassa virus; Lordsdale virus; Louping ill virus;Lymphocytic choriomeningitis virus; Machupo virus; Mayaro virus; MERScoronavirus; Measles virus; Mengo encephalomyocarditis virus; Merkelcell polyomavirus; Mokola virus; Molluscum contagiosum virus; Monkeypoxvirus; Mumps virus; Murray valley encephalitis virus; New York virus;Nipah virus; Norwalk virus; O'nyong-nyong virus; Orf virus; Oropouchevirus; Pichinde virus; Poliovirus; Punta toro phlebovirus; Puumalavirus; Rabies virus; Rift valley fever virus; Rosavirus A; Ross rivervirus; Rotavirus A; Rotavirus B; Rotavirus C; Rubella virus; Sagiyamavirus; Salivirus A; Sandfly fever sicilian virus; Sapporo virus; Semlikiforest virus; Seoul virus; Simian foamy virus; Simian virus 5; Sindbisvirus; Southampton virus; St. louis encephalitis virus; Tick-bornepowassan virus; Torque teno virus; Toscana virus; Uukuniemi virus;Vaccinia virus; Varicella-zoster virus; Variola virus; Venezuelan equineencephalitis virus; Vesicular stomatitis virus; Western equineencephalitis virus; WU polyomavirus; West Nile virus; Yaba monkey tumorvirus; Yaba-like disease virus; Yellow fever virus; Zika virus;

Anelloviridae. The Anelloviridae family consists of non-enveloped,circular, single-stranded DNA viruses. Three genera of anellovirus areknown to infect humans, named TTV, TTMDV, and TTMV.

Torque Teno Virus (TTV) is a non-enveloped, single-stranded DNA viruswith a circular, negative-sense genome. A smaller virus, which wassubsequently named Torque Teno-like Mini Virus (TTMV) has also beencharacterized, and a third virus with a genomic size in between that ofTTV and TTMV was discovered and subsequently named Torque Teno-like MidiVirus (TTMDV). Recent changes in nomenclature have classified the threeanelloviruses able to infect humans into Alphatorquevirus (TTV),Betatorquevirus (TTMV), and Gammatorquevirus (TTMDV) Genera of theAnelloviridae family of viruses. To date anelloviruses are stillconsidered “orphan” viruses waiting to be linked to human disease.

The human anelloviruses differ in genome size ranging from 3.8-3.9 kbfor TTV, 3.2 kb for TTMDV, and 2.8-2.9 kb for TTMV. A characteristicfeature of anelloviruses is the extreme diversity found both within andbetween anellovirus species; they can exhibit as much as 33%-50%divergence at the nucleotide level. Despite the nucleotide sequencediversity, anelloviruses share conserved genomic organization,transcriptional profiles, a non-coding GC rich region, and sequencemotifs resulting in shared virion structure and gene functions.

Anellovirus infections are highly prevalent in the general population. Astudy in Japan found that 75-100% of patients tested were infected withat least one of the three human anelloviruses, and many were infectedwith multiple species. Anelloviruses can infect young children, with theearliest documented infections occurring within the first months oflife. These viruses have been found in nearly every body site, fluid,and tissue tested including blood plasma, serum, peripheral bloodmononuclear cells (PBMCs), nasopharyngeal aspirates, bone marrow,saliva, breast milk, feces, as well as various tissues including thyroidgland, lymph node, lung, liver, spleen, pancreas, and kidney. Thereplication dynamics of anelloviruses are virtually unknown because ofthe inability to propagate these viruses in culture. Positive-strand TTVDNA, indicative of local viral replication, has been described inhepatocytes, bone marrow cells, and circulating PBMCs.

Anelloviruses are spread primarily through fecal-oral transmission,although mother-child and respiratory tract transmissions have also beenreported. There are conflicting reports regarding the presence of TTV incord blood specimens.

Reference sequences for anellovirus may be accessed at Genbank, e.g.Torque teno mini virus 1, Accession: NC_014097.1; Torque teno mini virus6, Accession: NC_014095.1; Torque teno midi virus 2, Accession:NC_014093.1; Torque teno midi virus 1, Accession: NC_009225.1; Torqueteno virus 3, Accession: NC_014081.1; Torque teno virus 19, Accession:NC_014078.1; Torque teno mini virus 8, Accession: NC_014068.1.

The term “antibiotic” as used herein includes all commonly usedbacteristatic and bactericidal antibiotics, usually those administeredorally. Antibiotics include aminoglycosides, such as amikacin,gentamicin, kanamycin, neomycin, streptomycin, and tobramycin;cephalosporins, such as cefamandole, cefazolin, cephalexin,cephaloglycin, cephaloridine, cephalothin, cephapirin, and cephradine;macrolides, such as erythromycin and troleandomycin; penicillins, suchas penicillin G, amoxicillin, ampicillin, carbenicillin, cloxacillin,dicloxacillin, methicillin, nafcillin, oxacillin, phenethicillin, andticarcillin; polypeptide antibiotics, such as bacitracin,colistimethate, colistin, polymyxin B; tetracyclines, such aschlortetracycline, demeclocycline, doxycycline, methacycline,minocycline, tetracycline, and oxytetracycline; and miscellaneousantibiotics such as chloramphenicol, clindamycin, cycloserine,lincomycin, rifampin, spectinomycin, vancomycin, and viomycin.Additional antibiotics'are described in “Remington's PharmaceuticalSciences,” 16th Ed., (Mack Pub. Co., 1980), pp. 1121-1178.

Antiviral agents. Individuals may receive antiviral therapy, which willalter the viral load for those viruses affected by the therapy. Examplesof viral infections thus treated include HIV, Bowenoid Papulosis,Chickenpox, Childhood HIV Disease, Human Cowpox, Hepatitis C, Dengue,Enteroviral, Epidermodysplasia Verruciformis, Erythema Infectiosum(Fifth Disease), Giant Condylomata Acuminata of Buschke and Lowenstein,Hand-Foot-and-Mouth Disease, Herpes Simplex, Herpes Virus 6, HerpesZoster, Kaposi Varicelliform Eruption, Rubeola Measles, Milker'sNodules, Molluscum Contagiosum, Monkeypox, Orf, Roseola Infantum,Rubella, Smallpox, Viral Hemorrhagic Fevers, Genital Warts, andNongenital Warts.

Antiviral agents include azidouridine, anasmycin, amantadine,bromovinyldeoxusidine, chlorovinyldeoxusidine, cytarbine, didanosine,deoxynojirimycin, dideoxycitidine, dideoxyinosine, dideoxynucleoside,desciclovir, deoxyacyclovir, edoxuidine, enviroxime, fiacitabine,foscamet, fialuridine, fluorothymidine, floxuridine, hypericin,interferon, interleukin, isethionate, nevirapine, pentamidine,ribavirin, rimantadine, stavirdine, sargramostin, suramin,trichosanthin, tribromothymidine, trichlorothymidine, vidarabine,zidoviridine, zalcitabine and 3-azido-3-deoxythymidine, and analogs,derivatives, pharmaceutically acceptable salts, esters, prodrugs,codrugs, and protected forms thereof.

Immunosuppression or immunosuppressive regimen, as used herein, refersto the treatment of an individual, for example a graft recipient withagents to diminish the immune responses of the host immune systemagainst autoantigens or graft. Exemplary immunosuppression regimens aredescribed in more detail herein.

Primary immunosuppressive agents include calcineurin inhibitors, whichcombine with binding proteins to inhibit calcineurin activity, and whichinclude, for example, tacrolimus, cyclosporine A, etc. Levels of bothcyclosporine and tacrolimus must be carefully monitored. Initially,levels can be kept in the range of 10-20 ng/mL, but, after 3 months,levels may be kept lower (5-10 ng/mL) to reduce the risk ofnephrotoxicity.

Adjuvant agents are usually combined with a calcineurin inhibitor andinclude steroids, azathioprine, mycophenolate mofetil, and sirolimus.Protocols of interest include a calcineurin inhibitor with mycophenolatemofetil. The use of adjuvant agents allows clinicians to achieveadequate immunosuppression while decreasing the dose and toxicity ofindividual agents. Mycophenolate mofetil in kidney transplant recipientshas assumed an important role in immunosuppression after severalclinical trials have shown a markedly decreased prevalence of acutecellular rejection compared with azathioprine and a reduction in 1-yeartreatment failures.

Antibody-based therapy may use monoclonal (eg, muromonab-CD3) orpolyclonal antibodies or anti-CD25 antibodies (eg, basiliximab,daclizumab) and is administered in the early post-transplant period (upto 8 wk). Antibody-based therapy allows for avoidance or dose reductionof calcineurin inhibitors, possibly reducing the risk of nephrotoxicity.The adverse effect profile of the polyclonal and monoclonal antibodieslimits their use in some patients.

The term “nucleic acid” used herein refers to a polynucleotidecomprising two or more nucleotides. It may be DNA or RNA. A “variant”nucleic acid is a polynucleotide having a nucleotide sequence identicalto that of its original nucleic acid except having at least onenucleotide modified, for example, deleted, inserted, or replaced,respectively. The variant may have a nucleotide sequence at least about80%, 90%, 95%, or 99%, identity to the nucleotide sequence of theoriginal nucleic acid.

Circulating, or cell-free, DNA was first detected in human blood plasmain 1948. (Mandel, P. Metais, P., C R Acad. Sci. Paris, 142, 241-243(1948)) Since then, its connection to disease has been established inseveral areas. (Tong, Y. K. Lo, Y. M., Clin Chim Acta, 363, 187-196(2006)) Studies reveal that much of the circulating nucleic acids inblood arise from necrotic or apoptotic cells (Giacona, M. B., et al.,Pancreas, 17, 89-97 (1998)) and greatly elevated levels of nucleic acidsfrom apoptosis is observed in diseases such as cancer. (Giacona, M. B.,et al., Pancreas, 17, 89-97 (1998); Fournie, G. J., et al., Cancer Lett,91, 221-227 (1995)). Particularly for cancer, where the circulating DNAbears hallmark signs of the disease including mutations in oncogenes,microsatellite alterations, and, for certain cancers, viral genomicsequences, DNA or RNA in plasma has become increasingly studied as apotential biomarker for disease. For example, Diehl et al recentlydemonstrated that a quantitative assay for low levels of circulatingtumor DNA in total circulating DNA could serve as a better marker fordetecting the relapse of colorectal cancer compared withcarcinoembryonic antigen, the standard biomarker used clinically.(Diehl, F., et al., Proc Natl Acad Sci, 102, 16368-16373 (2005); Diehl,F., et al., Nat Med, 14, 985-990 (2008)). Maheswaran et al reported theuse of genotyping of circulating cells in plasma to detect activatingmutations in epidermal growth factor receptors in lung cancer patientsthat would affect drug treatment. (Maheswaran, S., et al., N Engl J Med,359, 366-377 (2008)) These results collectively establish circulatingDNA free in plasma as a useful species in cancer detection andtreatment. Circulating DNA has also been useful in healthy patients forfetal diagnostics, with fetal DNA circulating in maternal blood servingas a marker for gender, rhesus D status, fetal aneuploidy, andsex-linked disorders. Fan et al recently demonstrated a strategy fordetecting fetal aneuploidy by shotgun sequencing of cell-free DNA takenfrom a maternal blood sample, a methodology that can replace moreinvasive and risky techniques such as amniocentesis or chorionic villussampling. (Fan, H. C., Blumenfeld, Y. J., Chitkara, U., Hudgins, L.,Quake, S. R., Proc Natl Acad Sci, 105, 16266-16271 (2008)).

The term “derived from” used herein refers to an origin or source, andmay include naturally occurring, recombinant, unpurified or purifiedmolecules. A nucleic acid derived from an original nucleic acid maycomprise the original nucleic acid, in part or in whole, and may be afragment or variant of the original nucleic acid. A nucleic acid derivedfrom a biological sample may be purified from that sample.

A “target nucleic acid” in the method according to the present inventionis a nucleic acid, DNA or RNA, to be detected. A target nucleic acidderived from an organism is a polynucleotide that has a sequence derivedfrom that of the organism and is specific to the organism. A targetnucleic acid derived from a pathogen refers to a polynucleotide having apolynucleotide sequence derived from that specific the pathogen.

In some embodiments, less than 1 pg, 5 pg, 10 pg, 20 pg, 30 pg, 40 pg,50 pg, 100 pg, 200 pg, 500 pg, 1 ng, 5 ng, 10 ng, 20 ng, 30 ng, 40 ng,50 ng, 100 ng, 200 ng, 500 ng, 1 μg, 5 μg, 10 μg, 20 μg, 30 μg, 40 μg,50 μg, 100 μg, 200 μg, 500 μg or 1 mg of nucleic acids are obtained fromthe sample for analysis. In some cases, about 1-5 pg, 5-10 pg, 10-100pg, 100 pg-1 ng, 1-5 ng, 5-10 ng, 10-100 ng, 100 ng-1 μg of nucleicacids are obtained from the sample for analysis.

In some embodiments, the methods described herein are used to detectand/or quantify nucleic acid sequences that correspond to a microbe ofinterest, or a microbiome of organisms. The methods described herein cananalyze at least 1; 2; 3; 4; 5; 10, 20; 50; 100; 200; 500; 1,000; 2,000;5,000; 10,000, 20,000; 50,000; 100,000; 200,000; 300,000; 400,000;500,000; 600,000; 700,000; 800,000; 900,000; 10⁶, 5×10⁶, 10⁷, 5×10⁷,10⁸, 5×10⁸, 10⁹ or more sequence reads.

In some embodiments, the methods described herein are used to detectand/or quantify gene expression, e.g. by determining the presence ofmRNA from a microorganism in relation to DNA from that microorganism. Insome embodiments, the methods described herein provide highdiscriminative and quantitative analysis of multiple genes. The methodsdescribed herein can discriminate and quantitate the expression of atleast 1, 2, 3, 4, 5, 10, 20, 50, 100, 200, 500, 1,000, 2,000, 5,000,10,000, 20,000, 50,000, 100,000, or more different target nucleic acids.

A sample containing cell-free nucleic acids is obtained from a subject.Such subject can be a human, a domesticated animal such as a cow,chicken, pig, horse, rabbit, dog, cat, goat, etc. In some embodiments,the cells used in the present invention are taken from a patient.Samples include, for example, the acellular fraction of whole blood,sweat, tears, saliva, ear flow, sputum, lymph, bone marrow suspension,lymph, urine, saliva, semen, vaginal flow, cerebrospinal fluid, brainfluid, ascites, milk, secretions of the respiratory, intestinal orgenitourinary tracts fluid, a lavage of a tissue or organ (e.g. lung) ortissue which has been removed from organs, such as breast, lung,intestine, skin, cervix, prostate, pancreas, heart, liver and stomach.Such samples can be separated by centrifugation, elutriation, densitygradient separation, apheresis, affinity selection, panning, FACS,centrifugation with Hypaque, etc. Once a sample is obtained, it can beused directly, frozen, or maintained in appropriate culture medium forshort periods of time.

To obtain a blood sample, any technique known in the art may be used,e.g. a syringe or other vacuum suction device. A blood sample can beoptionally pre-treated or processed prior to use. A sample, such as ablood sample, can be analyzed under any of the methods and systemsherein within 1 week, 6 days, 5 days, 4 days, 3 days, 2 days, 1 day, 12hrs, 6 hrs, 3 hrs, 2 hrs, or 1 hr from the time the sample is obtained.When obtaining a sample from a subject (e.g., blood sample), the amountcan vary depending upon subject size and the condition being screened.In some embodiments, up to 50, 40, 30, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2,or 1 mL of a sample is obtained. In some embodiments, 1-50, 2-40, 3-30,or 4-20 mL of sample is obtained. In some embodiments, more than 5, 10,15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95 or100 mL of a sample is obtained.

The cell-free fraction is preferably blood serum or blood plasma. Theterm “cell-free fraction” of a biological sample used herein refers to afraction of the biological sample that is substantially free of cells.The term “substantially free of cells” used herein refers to apreparation from the biological sample comprising fewer than about20,000 cells per ml, preferably fewer than about 2,000 cells per ml,more preferably fewer than about 200 cells per ml, most preferably fewerthan about 20 cells per ml. In contrast to certain prior art methods,genomic DNA is not excluded from the acellular sample, and typicallycomprises from about 50% to about 90% of the nucleic acids that arepresent in the sample.

The method of the present invention may further comprise preparing acell-free fraction from a biological sample. The cell-free fraction maybe prepared using conventional techniques known in the art. For example,a cell-free fraction of a blood sample may be obtained by centrifugingthe blood sample for about 3-30 min, preferably about 3-15 min, morepreferably about 3-10 min, most preferably about 3-5 min, at a low speedof about 200-20,000 g, preferably about 200-10,000 g, more preferablyabout 200-5,000 g, most preferably about 350-4,500 g. The biologicalsample may be obtained by ultrafiltration in order to separate the cellsand their fragments from a cell-free fraction comprising soluble DNA orRNA. Conventionally, ultrafiltration is carried out using a 0.22 μmmembrane filter.

The method of the present invention may further comprise concentrating(or enriching) the target nucleic acid in the cell-free fraction of thebiological sample. The target nucleic acid may be concentrated usingconventional techniques known in the art, such as solid phase absorptionin the presence of a high salt concentration, organic extraction byphenol-chloroform followed by precipitation with ethanol or isopropylalcohol, or direct precipitation in the presence of a high saltconcentration or 70-80% ethanol or isopropyl alcohol. The concentratedtarget nucleic acid may be at least about 2, 5, 10, 20 or 100 times moreconcentrated than that in the cell-free fraction. The target nucleicacid, whether or not concentrated, may be used for amplificationaccording to the method of the present invention.

In some embodiments the invention provides methods for diagnosis orprediction of transplant rejection. The term “transplant rejection”encompasses both acute and chronic transplant rejection. “Acuterejection or AR” is the rejection by the immune system of a tissuetransplant recipient when the transplanted tissue is immunologicallyforeign. Acute rejection is characterized by infiltration of thetransplanted tissue by immune cells of the recipient, which carry outtheir effector function and destroy the transplanted tissue. The onsetof acute rejection is rapid and generally occurs in humans within a fewweeks after transplant surgery. Generally, acute rejection can beinhibited or suppressed with immunosuppressive drugs such as rapamycin,cyclosporin A, anti-CD40L monoclonal antibody and the like.

“Chronic transplant rejection or CR” generally occurs in humans withinseveral months to years after engraftment, even in the presence ofsuccessful immunosuppression of acute rejection. Fibrosis is a commonfactor in chronic rejection of all types of organ transplants. Chronicrejection can typically be described by a range of specific disordersthat are characteristic of the particular organ. For example, in lungtransplants, such disorders include fibroproliferative destruction ofthe airway (bronchiolitis obliterans); in heart transplants ortransplants of cardiac tissue, such as valve replacements, suchdisorders include fibrotic atherosclerosis; in kidney transplants, suchdisorders include, obstructive nephropathy, nephrosclerorsis,tubulointerstitial nephropathy; and in liver transplants, such disordersinclude disappearing bile duct syndrome. Chronic rejection can also becharacterized by ischemic insult, denervation of the transplantedtissue, hyperlipidemia and hypertension associated withimmunosuppressive drugs.

In some embodiments, the invention further includes methods fordetermining the effectiveness of an immunosuppressive regimen for asubject who has received a transplant, e.g., an allograft.

Certain embodiments of the invention provide methods of predictingtransplant survival in a subject that has received a transplant. Theinvention provides methods of diagnosing or predicting whether atransplant in a transplant patient or subject will survive or be lost.In certain embodiments, the invention provides methods of diagnosing orpredicting the presence of long-term graft survival. By “long-term”graft survival is meant graft survival for at least about 5 years beyondcurrent sampling, despite the occurrence of one or more prior episodesof acute rejection. In certain embodiments, transplant survival isdetermined for patients in which at least one episode of acute rejectionhas occurred. As such, these embodiments provide methods of determiningor predicting transplant survival following acute rejection. Transplantsurvival is determined or predicted in certain embodiments in thecontext of transplant therapy, e.g., immunosuppressive therapy, whereimmunosuppressive therapies are known in the art. In yet otherembodiments, methods of determining the class and/or severity of acuterejection (and not just the presence thereof) are provided.

As in known in the transplantation field, the transplant organ, tissueor cell(s) may be allogeneic or xenogeneic, such that the grafts may beallografts or xenografts. A feature of the graft tolerant phenotypedetected or identified by the subject methods is that it is a phenotypewhich occurs without immunosuppressive therapy, i.e., it is present in ahost that is not undergoing immunosuppressive therapy such thatimmunosuppressive agents are not being administered to the host. Thetransplant graft maybe any solid organ and skin transplant. Examples oforgan transplants that can be analyzed by the methods described hereininclude but are not limited to kidney transplant, pancreas transplant,liver transplant, heart transplant, lung transplant, intestinetransplant, pancreas after kidney transplant, and simultaneouspancreas-kidney transplant.

Microbiome Detection and Analysis

The methods of the invention involve high throughput sequencing of acell-free nucleic acid sample from an individual, followed bybioinformatics analysis to determine the presence and prevalence ofmicrobial sequences, which sequences may be from indigenous organisms,e.g. the normal microbiome of gut, skin, etc., or may be non-indigenous,e.g. opportunistic, pathogenic, etc. infections. Analysis may beperformed for the complete microbiome, or for components there, forexample the virome, bacterial microbiome, fungal microbiome, protozoanmicrobiome, etc. Examples of nucleic acids include, but are not limitedto double-stranded DNA, single-stranded DNA, single-stranded DNAhairpins, DNA/RNA hybrids, RNA (e.g. mRNA or miRNA) and RNA hairpins. Insome embodiments, the nucleic acid is DNA. In some embodiments, thenucleic acid is RNA. For instance, cell-free RNA and DNA are present inhuman plasma.

Genotyping microbiome nucleic acids, and/or detection, identificationand/or quantitation of the microbiome-specific nucleic acids generallyinclude an initial step of amplification of the sample, although theremay be instances where sufficient cell free nucleic acids are availableand can be directly sequenced. When the nucleic acid is RNA, theamplification step may be preceded by a reverse transcriptase reactionto convert the RNA into DNA. Preferably the amplification is unbiased,that is the primers for amplification are universal primers, or adaptorsare ligated to the nucleic acids being analyzed, and amplificationprimers are specific for the adaptors. Examples of PCR techniquesinclude, but are not limited, to hot start PCR, nested PCR, in situpolonony PCR, in situ rolling circle amplification (RCA), bridge PCR,picotiter PCR and emulsion PCR. Other suitable amplification methodsinclude the ligase chain reaction (LCR), transcription amplification,self-sustained sequence replication, selective amplification of targetpolynucleotide sequences, consensus sequence primed polymerase chainreaction (CP-PCR), arbitrarily primed polymerase chain reaction(AP-PCR), degenerate oligonucleotide-primed PCR (DOP-PCR) and nucleicacid based sequence amplification (NABSA). Other amplification methodsthat may be used to amplify specific polymorphic loci include thosedescribed in, U.S. Pat. Nos. 5,242,794, 5,494,810, 4,988,617 and6,582,938.

Following amplification, the amplified nucleic acid is sequenced.Sequencing can be accomplished using high-throughput systems some ofwhich allow detection of a sequenced nucleotide immediately after orupon its incorporation into a growing strand, i.e., detection ofsequence in red time or substantially real time. In some cases, highthroughput sequencing generates at least 1,000, at least 5,000, at least10,000, at least 20,000, at least 30,000, at least 40,000, at least50,000, at least 100,000 or at least 500,000 sequence reads per hour;with each read being at least 50, at least 60, at least 70, at least 80,at least 90, at least 100, at least 120 or at least 150 bases per read.Sequencing can be performed using nucleic acids described herein such asgenomic DNA, cDNA derived from RNA transcripts or RNA as a template.

In some embodiments, high-throughput sequencing involves the use oftechnology available by Helicos BioSciences Corporation (Cambridge,Mass.) such as the Single

Molecule Sequencing by Synthesis (SMSS) method. SMSS is unique becauseit allows for sequencing an entire genome with no pre amplification stepneeded. Thus, distortion and nonlinearity in the measurement of nucleicacids are reduced. SMSS is described in part in US PublicationApplication Nos. 2006002471 I; 20060024678; 20060012793; 20060012784;and 20050100932.

In some embodiments, high-throughput sequencing involves the use oftechnology available by 454 Lifesciences, Inc. (Branford, Conn.) such asthe Pico Titer Plate device, which includes a fiber optic plate thattransmits chemiluminescent signal generated by the sequencing reactionto be recorded by a CCD camera in the instrument. This use of fiberoptics allows for the detection of a minimum of 20 million base pairs in4.5 hours.

Methods for using bead amplification followed by fiber optics detectionare described in Marguiles, M., et al. “Genome sequencing inmicrofabricated high-density pricolitre reactors”, Nature, doi:10.1038/nature03959; and well as in US Publication Application Nos.20020012930; 20030058629; 20030100102; 20030148344; 20040248161;20050079510, 20050124022; and 20060078909.

In some embodiments, high-throughput sequencing is performed usingClonal Single Molecule Array (Solexa, Inc.) or sequencing-by-synthesis(SBS) utilizing reversible terminator chemistry. These technologies aredescribed in part in U.S. Pat. Nos. 6,969,488; 6,897,023; 6,833,246;6,787,308; and US Publication Application Nos. 200401061 30;20030064398; 20030022207; and Constans, A., The Scientist 2003,17(13):36.

In some embodiments of this aspect, high-throughput sequencing of RNA orDNA can take place using AnyDot.chips (Genovoxx, Germany), which allowsfor the monitoring of biological processes (e.g., miRNA expression orallele variability (SNP detection). In particular, the AnyDot-chipsallow for 10×-50× enhancement of nucleotide fluorescence signaldetection. AnyDot.chips and methods for using them are described in partin International Publication Application Nos. WO 02088382, WO 03020968,WO 0303 1947, WO 2005044836, PCTEP 05105657, PCMEP 05105655; and GermanPatent Application Nos. DE 101 49 786, DE 102 14 395, DE 103 56 837, DE10 2004 009 704, DE 10 2004 025 696, DE 10 2004 025 746, DE 10 2004 025694, DE 10 2004 025 695, DE 10 2004 025 744, DE 10 2004 025 745, and DE10 2005 012 301.

Other high-throughput sequencing systems include those disclosed inVenter, J., et al. Science 16 Feb. 2001; Adams, M. et al, Science 24Mar. 2000; and M. J, Levene, et al. Science 299:682-686, January 2003;as well as US Publication Application No. 20030044781 and 2006/0078937.Overall such system involve sequencing a target nucleic acid moleculehaving a plurality of bases by the temporal addition of bases via apolymerization reaction that is measured on a molecule of nucleic acid,i.e., the activity of a nucleic acid polymerizing enzyme on the templatenucleic acid molecule to be sequenced is followed in real time. Sequencecan then be deduced by identifying which base is being incorporated intothe growing complementary strand of the target nucleic acid by thecatalytic activity of the nucleic acid polymerizing enzyme at each stepin the sequence of base additions. A polymerase on the target nucleicacid molecule complex is provided in a position suitable to move alongthe target nucleic acid molecule and extend the oligonucleotide primerat an active site. A plurality of labeled types of nucleotide analogsare provided proximate to the active site, with each distinguishablytype of nucleotide analog being complementary to a different nucleotidein the target nucleic acid sequence. The growing nucleic acid strand isextended by using the polymerase to add a nucleotide analog to thenucleic acid strand at the active site, where the nucleotide analogbeing added is complementary to the nucleotide of the target nucleicacid at the active site. The nucleotide analog added to theoligonucleotide primer as a result of the polymerizing step isidentified. The steps of providing labeled nucleotide analogs,polymerizing the growing nucleic acid strand, and identifying the addednucleotide analog are repeated so that the nucleic acid strand isfurther extended and the sequence of the target nucleic acid isdetermined.

In some embodiments, shotgun sequencing is performed. In shotgunsequencing, DNA is broken up randomly into numerous small segments,which are sequenced using the chain termination method to obtain reads.Multiple overlapping reads for the target DNA are obtained by performingseveral rounds of this fragmentation and sequencing. Computer programsthen use the overlapping ends of different reads to assemble them into acontinuous sequence.

In some embodiments, the invention provides methods for detection andquantitation of microbial sequences using sequencing. In this case, onecan estimate the sensitivity of detection. There are two components tosensitivity: (i) the number of molecules analyzed (depth of sequencing)and (ii) the error rate of the sequencing process. Regarding the depthof sequencing, a frequent estimate for the variation between individualsis that about one base per thousand differs. Currently, sequencers suchas the Illumina Genome Analyzer have read lengths exceeding 36 basepairs. While the fraction of host DNA in the blood may be variabledepending on the status of the individual, one can take 90% as abaseline estimate. At this fraction of donor DNA, approximately one in10 molecules analyzed will be microbial. On the Genome Analyzer one canobtain about 10 million molecules per analysis channel and there are 8analysis channels per instrument run. Therefore, if one sample is loadedper channel, one should be able to detect about 10⁶ molecules that canbe identified as microbial and informative as to the status of themicrobiome. Higher sensitivity can be achieved simply by sequencing moremolecules, i.e. using more channels.

The sequencing error rate also affects the sensitivity of thistechnique. Typical sequencing error rates for base substitutions varybetween platforms, but are between 0.5-1.5%. This places a potentiallimit on sensitivity of 0.16 to 0.50%. However, it is possible tosystematically lower the sequencing error rate by resequencing thesample template multiple times, as has been demonstrated by HelicosBioSciences (Harris, T. D., et al., Science, 320, 106-109 (2008)). Asingle application of resequencing would reduce the expected error rate.

Following sequencing, the dataset of sequences is uploaded to a dataprocessor for bioinformatics analysis to subtract host sequences, i.e.human, cat, dog, etc. from the analysis; and determine the presence andprevalence of microbial sequences, for example by a comparison of thecoverage of sequences mapping to a microbial reference sequence tocoverage of the host reference sequence. The subtraction of hostsequences may include the step of identifying a reference host sequence,and masking microbial sequences or microbial-mimicking sequences presentin the reference host genome. Similarly, determining the presence of amicrobial sequence by comparison to a microbial reference sequence mayinclude the step of identifying a reference microbial sequence, andmasking host sequences or host-mimicking sequences present in thereference microbial genome.

The dataset is optionally cleaned to check sequence quality, removeremnants of sequencer specific nucleotides (adapter sequences), andmerge paired end reads that overlap to create a higher quality consensussequence with less read errors. Repetitive sequences are identified asthose having identical start sites and length, and duplicates may beremoved from the analysis.

An important feature of the invention is the subtraction of humansequences from the analysis. As the amplification/sequencing steps areunbiased, the preponderance of sequences in a sample will be hostsequences. The subtraction process may be optimized in several ways toimprove the speed and accuracy of the process, for example by performingmultiple subtractions where the initial alignment is set at a coarsefilter, i.e. with a fast aligner, and performing additional alignmentswith a fine filter, i.e. a sensitive aligner.

The database of reads are initially aligned against a human referencegenome, including without limitation Genbank hg19 reference sequences,to bioinformatically subtract the host DNA. Each sequence is alignedwith the best fit sequence in the human reference sequence Sequencespositively identified as human are bioinformatically removed from theanalysis.

The reference human sequence can also be optimized by adding in contigsthat have a high hit rate, including without limitation highlyrepetitive sequence present in the genome that are not well representedin reference databases. It has been observed that of the reads that donot align to hg19, a significant amount is eventually identified ashuman in a later stage of the pipeline, when a more sensitive aligner isused. Removing these reads earlier in the analysis can be performed bybuilding an expanded human reference. This reference is created byidentifying human contigs in a human sequence database other than thereference, e.g. NCBI NT database, that have high coverage after theinitial human read subtraction. Those contigs are added to the humanreference to create a more comprehensive reference set. Additionallynovel assembled human contigs from cohorts studies can be used as afurther mask for human-derived reads.

Regions of the human genome reference sequence that contain non-humansequences may be masked, e.g. viral and bacterial sequences that areintegrated into the genome of the reference sample. For example,Epstein-Barr Virus (EBV) has about 80% of its genome incorporated intohg19.

Sequence reads identified as non-human are then aligned to a nucleotidedatabase of microbial reference sequences. The database may be selectedfor those microbial sequences known to be associated with the host, e.g.the set of human commensal and pathogenic microorganisms.

The microbial database may be optimized to mask or remove contaminatingsequences. For example it has been observed that many public databaseentries include artifactual sequences not derived from themicroorganism, e.g., primer sequences, host sequences, and othercontaminants. It is desirable to perform an initial alignment orplurality of alignments on a database. Regions that show irregularitiesin read coverage when multiple samples are aligned can be masked orremoved as an artifact. The detection of such irregular coverage can bedone by various metrics, such as the ratio between coverage of aspecific nucleotide and the average coverage of the entire contig withinwhich this nucleotide is found. In general a sequence that isrepresented as greater than about 5×, about 10×, about 25×, about 50×,about 100× the average coverage of that reference sequence areartifactual. Alternatively a binomial test can be applied to provide aper-base likelihood of coverage given the overall coverage of thecontig. Removal of contaminant sequence from reference databases allowsaccurate identification of microbes. It is a benefit of the methods ofthe invention that the databases are improved with alignment of samples,e.g. a database my be aligned with 1, 10, 20, 50, 100 or more samples toimprove the database prior to commercial or clinical use.

Each high confidence read may align to multiple organisms in the givenmicrobial database. To correctly assign organism abundance based uponthis possible mapping redundancy, an algorithm is used to compute themost likely organism an algorithm is selected (for example see Lindneret al. Nucl. Acids Res. (2013) 41 (1): e10). For example GRAMMy or GASiCalgorithms can be used to compute the most likely organism that a givenread came from. These data provide information regarding the presence ofa microbe in the cell-free nucleic acid sample.

These assignments are then totaled and used to compute the estimatednumber of reads assigned to each organism in a given sample in adetermination of the prevalence of the organism in the cell-free nucleicacid sample. The analysis normalizes the counts for the size of themicrobial genome to provide a calculation of coverage for the microbe.The normalized coverage for each microbe is compared to the hostsequence coverage in the same sample to account for differences insequencing depth between samples.

The final determination provides a dataset of microbial organismsrepresented by sequences in the sample, and the prevalence of thosemicroorganisms. These data are optionally aggregated and displayed forready visualization, e.g. in the form of a report provided to theindividual or health care provider; or written in a browser format withhyper-linked data. The coverage estimation can be aggregated withmetadata from the sample, and sorted into tables and figures for eachsample, or cohort of samples.

Optionally the host sequences that are filtered out can be used forother purposes, e.g. in personalized medicine. For example, certain SNPsin the human genome may allow doctors to identify drug sensitivities fora given patient. The human-derived sequences may reveal integration ofviruses into the host's genome (e.g. EBV, HPV, poliomavirus) or be usedfor synergistic clinical applications (e.g., cell-free tumor DNA may beused to monitor cancer progressing in parallel with infection monitoringin patients that are highly susceptible to infection due tochemotherapy).

Also provided are reagents and kits thereof for practicing one or moreof the above-described methods. The subject reagents and kits thereofmay vary greatly. Reagents of interest include reagents specificallydesigned for use in production of the above-described: (i) profiling ofa microbiome and an individual; (ii) identification of microbiomeprofiles; and (ii) detection and/or quantitation of one or more nucleicacids from a microbiome in a sample obtained from an individual. Thekits may comprise reagents necessary to perform nucleic acid extractionand/or nucleic acid detection using the methods described herein such asPCR and sequencing. The kit may further comprise a software package fordata analysis, which may include reference profiles for comparison withthe test profile., and in particular may include reference databasesoptimized as described above. The kits may comprise reagents such asbuffers, and H₂O.

Such kits may also include information, such as scientific literaturereferences, package insert materials, clinical trial results, and/orsummaries of these and the like, which indicate or establish theactivities and/or advantages of the composition, and/or which describedosing, administration, side effects, drug interactions, or otherinformation useful to the health care provider. Such kits may alsoinclude instructions to access a database. Such information may be basedon the results of various studies, for example, studies usingexperimental animals involving in vivo models and studies based on humanclinical trials. Kits described herein can be provided, marketed and/orpromoted to health providers, including physicians, nurses, pharmacists,formulary officials, and the like. Kits may also, in some embodiments,be marketed directly to the consumer.

Any of the methods above can be performed by a computer program productthat comprises a computer executable logic that is recorded on acomputer readable medium. For example, the computer program can executesome or all of the following functions: (i) controlling isolation ofnucleic acids from a sample, (ii) pre-amplifying nucleic acids from thesample, (iii) amplifying, sequencing or arraying specific regions in thesample, (iv) identifying and quantifying a microbial sequence in thesample, (v) comparing data on a microbe presence or prevalence detectedfrom the sample with a predetermined threshold, (vi) determininginfection, microbiome health, immunocompetence status or outcome, (vi)declaring the sample status with respect to infection, microbiomehealth, immunocompetence, etc.

The computer executable logic can work in any computer that may be anyof a variety of types of general-purpose computers such as a personalcomputer, network server, workstation, or other computer platform now orlater developed. In some embodiments, a computer program product isdescribed comprising a computer usable medium having the computerexecutable logic (computer software program, including program code)stored therein. The computer executable logic can be executed by aprocessor, causing the processor to perform functions described herein.In other embodiments, some functions are implemented primarily inhardware using, for example, a hardware state machine. Implementation ofthe hardware state machine so as to perform the functions describedherein will be apparent to those skilled in the relevant arts.

The program can provide a method of evaluating microbial status in anindividual by accessing data that reflects the profiling of themicrobiome and the individual, and/or the quantitation of one or morenucleic acids from the microbiome in the circulation of the individual.

In one embodiment, the computer executing the computer logic of theinvention may also include a digital input device such as a scanner. Thedigital input device can provide information on a nucleic acid, e.g.,presence or prevalence

In some embodiments, the invention provides a computer readable mediumcomprising a set of instructions recorded thereon to cause a computer toperform the steps of (i) receiving data from one or more nucleic acidsdetected in a sample; and (ii) diagnosing or predicting a status basedon the microbiome quantitation.

Also provided are databases of microbial reference sequences; anddatabases of human reference sequences. Such databases will typicallycomprise optimized datasets as described above.

In some embodiments the methods of the invention provide an individual'sstatus with respect to infection. In some such embodiments, themicrobial infection is a pathogen, where any presence of the pathogensequence indicates a clinically relevant infection. In otherembodiments, the prevalence is indicative of microbial load, where apre-set level is indicative of clinical relevance. In some suchembodiments the individual is treated or considered for treatment withan antimicrobial therapy, e.g. antibiotics, passive or activeimmunotherapy, antivirals, etc. An individual may be tested beforetherapy, during therapy, and after therapy.

A microbial infection may also be indicated by load for a commensalorganism, where the level of a commensal in a blood sample is indicativeof intestinal health, e.g. gut lumen breakdown.

A comparison can be made of microbioal RNA, alone or in relation tomicrobiol DNA, where an excess of RNA for a microbial sequence, e.g.about 5×, 10×, 15×, 20×, 25× the coverage of the microbial DNA isindicative of an active infection. In some embodiments the microbe thusanalyzed is one capable of latent infection, e.g. herpesvirus, hepatitisvirus, etc.

In other embodiments, an overall estimate of the microbiome is ofinterest, where the relative presence of prevalence of classes ofmicroorganisms are of interest. It is known in the art that diet andtreatment with drugs, e.g. statins, antibiotics, immunosuppressiveagents, etc. can affect the overall health of the microbiome, and it istherefor of interest to determine the composition of the microbiome.

In some embodiments, temporal differences in the amount of said one ormore nucleic acids from the microbiome can be used to monitoreffectiveness of anti-microbial treatment or to select an treatment. Forinstance, the amount of one or more nucleic acids from the microbiomecan be determined before and after an treatment. A decrease in the oneor more nucleic acids from the microbe after treatment may indicate thatthe treatment was successful. Additionally, the amount of one or morenucleic acids from the microbiome can be used to choose betweentreatments, for examples, treatments of different strengths.

In one aspect the invention provides methods for the diagnosis orprediction of immunocompetence, transplant status or outcome in asubject receiving an immunosuppressive regimen. Followingimmunosuppression, samples as described above can be drawn from thepatient and analyzed for the presence or absence of one or moremicrobiome, including virome nucleic acids. In some embodiments, thesample is blood, plasma, serum or urine. The proportion and/or amount ofmicrobial nucleic acids can be monitored over time and an increase inthis proportion can be used to determine immunocompetence. Thequantitation of load may be determined by any suitable method known inthe art including those described herein such as sequencing, nucleicacid arrays or PCR.

In some embodiments, the amount of one or more microbiome nucleic acidsin a sample from the immunosuppressed recipient is used to determine thetransplant status or outcome. Thus, in some embodiments, the methods ofthe invention further comprise quantitating the one or more nucleicacids from the microbiome. In some embodiments, the amount of one ormore nucleic acids from the donor sample is determined as a percentageof total the nucleic acids in the sample. In some embodiments, theamount of one or more nucleic acids from the donor sample is determinedas a ratio of the total nucleic acids in the sample. In someembodiments, the amount of one or more nucleic acids from the donorsample is determined as a ratio or percentage compared to one or morereference nucleic acids in the sample. For instance, the amount of oneor more nucleic acids from the microbiome can be determined to be 10% ofthe total nucleic acids in the sample. Alternatively, the amount of oneor more nucleic acids from the microbiome can be at a ratio of 1:10compared to total nucleic acids in the sample. Further, the amount ofone or more nucleic acids from the microbiome can be determined to be10% or at a ratio of 1:10 of a reference gene such a β-globin. In someembodiments, the amount of one or more nucleic acids from the microbiomecan be determined as a concentration. For example, the amount of one ormore nucleic acids from the donor sample can be determined to be 1μg/mL.

In some embodiments, the amount of one or more nucleic acids from themicrobiome above a predetermined threshold value is indicative of aimmunocompetence status. For example, the normative values forclinically stable patients with no evidence of graft rejection or otherpathologies can be determined. An increase in the amount of one or morenucleic acids from the microbiome below the normative values forclinically stable post-transplantation patients could indicate a stableoutcome. On the other hand, an amount of one or more nucleic acids fromthe microbiome above or at the normative values for clinically stablepost-transplantation patients could indicate increased immunocompetenceand risk of graft rejection.

In some embodiments, different predetermined threshold values areindicative of different transplant outcomes or status. For example, asdiscussed above, an increase in the amount of one or more nucleic acidsfrom the microbiome above the normative values for clinically stablepost-transplantation patients could indicate a change in transplantstatus or outcome such as transplant rejection or transplant injury.However, an increase in the amount of one or more nucleic acids from themicrobiome above the normative values for clinically stablepost-transplantation patients but below a predetermined threshold levelcould indicate a less serious condition such as a viral infection ratherthan transplant rejection. An increase in the amount of one or morenucleic acids from the microbiome above a higher threshold couldindicate transplant rejection.

In some embodiments, temporal differences in the amount of said one ormore nucleic acids from the microbiome are indicative ofimmunocompetence. For instance, a transplant patient can be monitoredover time to determine the amount of one or more nucleic acids from themicrobiome. A temporary decrease in the amount of one or more nucleicacids from the microbiome, which subsequently return to normal values,might indicate a less serious condition rather than transplantrejection. On the other hand, a sustained decrease in the amount one ormore nucleic acids from the microbiome might indicate a seriouscondition such as lack of effective immunosuppression and graftrejection.

In some embodiments, temporal differences in the amount of said one ormore nucleic acids from the microbiome can be used to monitoreffectiveness of an immunosuppressant treatment or to select animmunosuppressant treatment. For instance, the amount of one or morenucleic acids from the microbiome can be determined before and after animmunosuppressant treatment. A decrease in the one or more nucleic acidsfrom the microbiome after treatment may indicate that the treatment wassuccessful in preventing transplant rejection. Additionally, the amountof one or more nucleic acids from the microbiome can be used to choosebetween immunosuppressant treatments, for examples, immunosuppressanttreatments of different strengths. For example, a lower amount in one ormore nucleic acids from the microbiome may indicate that there is a needof a very potent immunosuppressant, whereas a higher amount in one ormore nucleic acids from the microbiome may indicate that a less potentimmunosuppressant may be used.

The invention provides methods that sensitive and specific. In someembodiments, the methods described herein for diagnosing or predictingtransplant status or outcome have at least 56%, 60%, 70%, 80%, 90%, 95%or 100% sensitivity. In some embodiments, the methods described hereinhave at least 56% sensitivity. In some embodiments, the methodsdescribed herein have at least 78% sensitivity. In some embodiments, themethods described herein have a specificity of about 70% to about 100%.In some embodiments, the methods described herein have a specificity ofabout 80% to about 100%. In some embodiments, the methods describedherein have a specificity of about 90% to about 100%. In someembodiments, the methods described herein have a specificity of about100%.

The invention provides non-invasive diagnostics for individuals,including individuals that are being treated with immunosuppressiveregimens, treated with anti-microbial agents, etc., by monitoring thesequences of cell-free DNA or RNA from non-human sources. For example,individuals carry a number of virus, where the virus load is shownherein to vary with the immunocompetence of the individual. Preferredvirus for monitoring immunocompetence are annellovirus, in which theviral burden is shown herein to correlate with immunocompetence of theindividual.

In some embodiments, the invention provides methods, devices,compositions and kits for detection and/or quantitating circulatingnucleic acids, usually free in plasma or from viral particles, for thediagnosis, prognosis, detection and/or treatment of an infection, ofimmunocompetence, transplant status or outcome.

In some specific embodiments, the invention provides an approach tononinvasive detection of immunocompetence in transplant patients byvirome analysis, which circumvents the potential problems ofmicrochimerism from DNA from other foreign sources and is general forall organ recipients without consideration of gender. In someembodiments, a genetic fingerprint is generated for the virome of theindividual. This approach allows for a reliable identification ofsequences that can be made in a manner that is independent of thegenders of donor and recipient.

Following an immunosuppressive regimen, e.g. in conjunction withtransplantation, treatment of autoimmune disease, etc., bodily fluidsuch as blood can be drawn from the patient and analyzed for markers.Examples of bodily fluids include, but are not limited to, smears,sputum, biopsies, secretions, cerebrospinal fluid, bile, blood, lymphfluid, saliva, and urine. Detection, identification and/or quantitationof the virome sequences can be performed using real-time PCR, chips,high-throughput shotgun sequencing of circulating nucleic acids (e.g.cell-free DNA), as well as other methods known in the art including themethods described herein. The viral load can be monitored over time andan increase in this proportion can be used to determine immunocompetencestatus or outcome.

In any of the embodiments described herein, the transplant graft can beany solid organ or skin transplant. Examples of organ transplants thatcan be analyzed by the methods described herein include but are notlimited to kidney transplant, pancreas transplant, liver transplant,heart transplant, lung transplant, intestine transplant, pancreas afterkidney transplant, and simultaneous pancreas-kidney transplant.

In some other embodiments, the methods of the invention are used indetermining the efficacy of a therapy for treatment of disease,including infection, either at an individual level, or in the analysisof a group of patients, e.g. in a clinical trial format. Suchembodiments typically involve the comparison of two time points for apatient or group of patients. The patient status is expected to differbetween the two time points as the result of a therapeutic agent,therapeutic regimen, or disease challenge to a patient undergoingtreatment.

Examples of formats for such embodiments may include, withoutlimitation, analyzing the microbiome at two or more time points, where afirst time point is a diagnosed but untreated patient; and a second oradditional time point(s) is a patient treated with a candidatetherapeutic agent or regimen.

In another format, a first time point is a diagnosed patient in diseaseremission, e.g. as ascertained by current clinical criteria, as a resultof a candidate therapeutic agent or regimen. A second or additional timepoint(s) is a patient treated with a candidate therapeutic agent orregimen, and challenged with a disease-inducing agent, for example in avaccine context.

In such clinical trial formats, each set of time points may correspondto a single patient, to a patient group, e.g. a cohort group, or to amixture of individual and group data. Additional control data may alsobe included in such clinical trial formats, e.g. a placebo group, adisease-free group, and the like, as are known in the art. Formats ofinterest include crossover studies, randomized, double-blind,placebo-controlled, parallel group trial is also capable of testing drugefficacy, and the like. See, for example, Clinical Trials: AMethodologic Perspective Second Edition, S. Piantadosi,Wiley-Interscience; 2005, ISBN-13: 978-0471727811; and Design andAnalysis of Clinical Trials: Concepts and Methodologies, S. Chow and J.Liu, Wiley-Interscience; 2003; ISBN-13: 978-0471249856, each hereinspecifically incorporated by reference.

EXAMPLES Temporal Response of the Human Virome to Immunosuppression andAntiviral Therapy

The viral component of the microbiome, the human virome, remainsrelatively understudied (Wylie et al. (2012) Transl Res 160, 283-290)and little is known about the effects of immune modulation and antiviraltherapies on virome composition. It was previously shown that thehealthy gut virome remains remarkably stable over time (Reyes et al.(2010) Nature 466, 334-338), and that the predominant source ofvariation is due to differences between subjects, although anassociation between diet and the virome composition was found (Minot, etal. (2011). Genome Research 21, 1616-1625).

Immunosuppressive therapies significantly reduce the risk of graftrejection in organ transplantation but increase the susceptibility ofrecipients to infections. Infections with viral pathogens, in particularthe herpesvirus cytomegalovirus (CMV), occur frequently and increase therecipient's risk of graft failure. Organ transplant recipients aretherefore frequently subjected to antiviral prophylactic or preemptivetherapies directed against CMV.

The inverse relationship between the level of immunosuppression and therisks of infection and rejection leaves only a narrow therapeutic windowavailable for patient treatment. Post-transplant care is furthercomplicated by numerous limitations of the currently available methodsfor the diagnosis of infection and rejection. Diagnosis of rejectionmostly relies on invasive biopsies that suffer from interobservervariability, high cost and patient discomfort. Diagnosis of infectionsis challenging given the fact that the symptoms of infection arediminished following immunosuppression, and commonly used diagnosticmethods, such as antigen-detection and PCR-based molecular tests, relyon a specific target and therefore an a priori hypothesis for the sourceof the infection.

As a final complication, patient-to-patient variability in thesensitivity to immunosuppressive drugs can give rise to over- and underimmunosuppression, increasing the risk of infection or rejectionrespectively.

There are few substantive methods to measure the health of the immunesystem, and the connection between immunocompetence and the viralcomponent of the microbiome is poorly understood. Organ transplantrecipients are treated with a post-transplant therapy that combinesimmunosuppressive and antiviral drugs, offering a window into theeffects of immune modulation on the human virome. We used sequencing ofcell-free DNA in plasma to investigate drug-virome interactions in acohort of organ transplant recipients (656 samples, 96 patients), andfind that antivirals and immunosuppressants strongly affect thestructure of the virome in plasma. We observe marked viromecompositional dynamics at the onset of the therapy and find that thetotal viral load increases with immunosuppression, whereas the bacterialcomponent of the microbiome remains largely unaffected. The data provideinsight into the relationship between the human virome, the state of theimmune system, and the effects of pharmacological treatment, and offer apotential application of the virome state to predict immunocompetence.

In this work, we sequenced cell-free DNA circulating in plasma toinvestigate drug-microbiome interactions following organtransplantation. We studied the patterns of infection in heart and lungtransplant recipients subjected to a combination of immunosuppressantsand antiviral prophylaxis. We find that immunosuppressants andantivirals have a strong influence on the structure of the viralcomponent of the microbiome but not the bacterial component. Strongcompositional dynamics are observed at the onset of the drug therapy asthe virome composition of different individuals converge to a similar,drug-determined state. The total viral load increases markedly inresponse to the therapy, as viruses, in particular the anelloviruses,take advantage of a reduction of immunocompetence. Finally, we show thatmeasurement of the anellovirus burden enables stratification ofrejecting and non-rejecting recipients.

656 plasma samples were collected longitudinally from 96 solid organtransplant recipients (41 adult heart, 24 pediatric heart, 31 adultlung). Cell-free DNA was purified from plasma and sequenced. In total,we obtained 820 gigabases (Gbp) of sequencing data, with an average of1.25 Gbp per sample (Illumina HiSeq, 1×50 bp reads, FIG. 1B). Organtransplant recipients were continuously enrolled in the study over thecourse of more than 2 years and samples were collected from therecipients at regular time points post transplant, with the highestfrequency of sample collection in the first months post transplant. FIG.1C shows the number of samples analyzed as a function of time posttransplant for the different patient classes.

The patients in the cohort were treated with antiviral prophylaxis andimmunosuppression as part of a standardized post-transplant therapy(FIG. 1D). Maintenance immunosuppression was tacrolimus-based for theadult heart and lung transplant recipients and was complemented withmycophenolate mofetil and prednisone. Pediatric patients were treatedwith a cyclosporine based anti-rejection therapy. CMV positivetransplant recipients (prior CMV infection for recipient and/or donor),but not CMV negative recipients, were treated with antiviralprophylaxis. The protocol design entails high doses ofimmunosuppressants and antiviral drugs in the first few months posttransplant, after which the doses are gradually reduced as the risks ofrejection and infection diminish. Given the narrow therapeutic windowavailable for immunosuppression and the large patient-to-patientvariability in pharmacokinetics of tacrolimus, the concentration of thetacrolimus is directly measured in the blood and the dose is adjusted tomaintain a target drug level. FIG. 1D shows the mean level of tacrolimusmeasured in blood for the tacrolimus-treated patients and illustratesthe design of the drug treatment protocol.

DNA sequence analysis. Microbiome-derived sequences were identifiedafter computational subtraction of human-derived sequences. To this end,duplicate and low quality reads were removed and the remaining readswere mapped to the human reference genome, build hg19 (BWA (Li andDurbin, 2009), see methods). Unmapped reads were then collected and lowcomplexity reads were removed. FIG. 1E shows the distribution of theremaining read fraction after applying duplicate and quality filters(average of 86%) and the distribution of the remaining fraction aftersubtraction of human reads (average of 2%).

To identify infectious agents, the remaining, high quality, unique,non-human reads were mapped using BLAST to a reference database of viral(n=1401), bacterial (n=1980) and fungal (n=32) genomes (downloaded fromNCBI, FIG. 6A). 0.12% of the uniquely sequenced reads aligned to atleast one of the target genomes (FIG. 6B, C). We used a quantitative PCR(qPCR) assay targeted to a subset of sequencing identified targets(herpesviruses 4,5,6 and parvovirus) to validate the positive hitsidentified by the sequencing-based approach. We found a quantitativeagreement between viral counts as measured by sequencing and qPCR (FIG.6D).

We furthermore found that the sensitivity of the sequencing assay forthe detection of herpesviruses is on par with qPCR measurements. Thelarger capture cross-section available to the sequencing assay—thecomplete target genome versus the PCR amplicon target region—is thussufficient to overcome the signal loss in sequencing, caused by thefinite efficiency of sequencing library preparation and libraryundersampling. The highest CMV loads measured using sequencing acrossall samples in the study corresponded to two adult heart transplantpatients that suffered from a clinically diagnosed disseminated CMVinfection (see FIG. 6E).

To test for the presence of potential contaminants in the reagents usedfor DNA extraction and sequencing library preparation, we performed twocontrol experiments. In the first, we prepared 2 samples with a knowntemplate (Lambda gDNA, Pacbio Part no: 001-119-535), and purified DNAfor sequencing using the above-described workflow (Illumina Miseq, 3.4and 3.5 million reads). Lambda-derived sequences were removed and theremaining sequences (0.4%) were aligned to the BLAST reference databasedescribed above. No evidence was found for the various infectious agentsdiscussed in this work, but we did detect sequences related to theEnterobacteriaceae bacterial family (phylum Proteobacteria), primarilyE. coli (>97%), and enterobacterial phages (<1%), which are likely aremnant of the lambda DNA culture. In a second control, we prepared asample for sequencing from nuclease-free water. The sample was includedin a sequencing run along with a sample unrelated to this work andrecruited only a limited number of sequences, 15 in total, which mappedto genomes of two bacterial species. Again, no evidence was found forthe infectious agents that are discussed below.

We studied the microbiome composition in plasma at different levels oftaxonomic classification using Grammy, a tool that utilizes thesequence-similarity data obtained with BLAST to perform a maximumlikelihood estimation of the relative abundance of species. GRAMMyaccounts for differences in target genome size and the ambiguity of readassignments. Note that this approach only allows estimating theabundance of species for which genomic data is available in thereference database. FIG. 1F shows the relative abundance of species atdifferent levels of taxonomic classification (average over all samples).We find that viruses (73%) are more abundantly represented than bacteria(25%) and fungi (2%) (FIG. 1F panel a). Among viruses, we find thatssDNA viruses occupy a larger fraction (72%) than dsDNA viruses (28%).Seven distinct viral families are found (abundance>0.75%), with onedominant family, the Anelloviridae, which accounted for 68% of the totalpopulation (FIG. 1F, panel b). The anelloviridae fraction is mostly(97%) composed of viruses from the Alphatorquevirus genus (FIG. 1F,panel c). The Alphatorque genus is the genus of Torque Teno Viruses(TTVs), and sequences related to 14 different torque teno virotypes wereidentified (FIG. 1, panel d). Infections with polyomaviruses arewidespread in the human population, and polyomavirus DNAemia is notuncommon in the first year after solid organ transplantation.Polyomavirus-derived sequences were found in 75 samples (11%)corresponding to 36 patients in the present cohort. Evidence for thepresence of BK (41%), JC (27%), TS (4%), WU polyomavirus (6%), SV40 (6%)and the recently discovered HPyV6 (13%) (Schowalter et al., 2010) wasfound (FIG. 1F, panel e). Among bacteria, Proteobacteria (36%),Firmicutes (50%), Actinobacteria (10%), Bacteroidetes (4%) are the phylamost abundantly represented in the sample (FIG. 1F, panel f).

To investigate potential incorrect assignments of the relatively shortreads available to this study (50 bp), we examined the dependence of theabundance estimates on read length, based on longer, paired-end reads(2×100 bp) collected for a subset of samples (n=55). We found that theabundance estimates based on 50 bp subreads and 100 bp reads are similarfor all levels of taxonomic classification reported here (FIG. 6F).

Sensitivity of virome composition to drug dosage. The available clinicaldata on drug dosage was used to analyze drug-microbiome interactions.Here, we examined data for the adult heart and lung transplant patientsthat were treated with a tacrolimus-based anti-rejection protocol (47patients and 380 observations), thereby excluding the pediatric patientsthat were treated with cyclosporine and patients that were switched fromtacrolimus to cyclosporine immunosuppression due to drug-intoleranceissues. Data on prescription antiviral drug doses (valganciclovir) andthe measured levels of tacrolimus in blood were collected fromindividual patient records and the mean composition for samplescorresponding to different drugs levels was extracted. To account for adelayed effect of the microbiome composition on dose changes, the druglevel and dose data were sliding window average filtered (see FIG. 1Cand FIG. 7A-C; window size 45 days).

We find that the structure of the viral component of the microbiome is asensitive function of drug dosage (47 patients, 380 samples, FIG. 2A).However, the structure of the bacterial component of the microbiome wasnot significantly altered by the drug therapy, as discussed furtherbelow (FIG. 7D). Herpesvirales and caudovirales dominated the viromewhen patients received a low dose of valganciclovir and tacrolimus. Incontrast, a high dose of immunosuppressants and antivirals gave rise toa virome structure that is dominated by anelloviridae (up to 94%occupation at high drug levels). The antiviral prophylaxis is intendedto prevent CMV disease, but other herpesviruses are also susceptible tothe drug so it is not surprising that a higher dose of valganciclovirgives rise to a lower fraction of viruses from the Herpesvirales order.The observation that anelloviridae take advantage of suppression of thehost immune system is consistent with various observations from theliterature: it was previously shown that the incidence of anelloviridaeincreases with progression towards AIDS in HIV patients, and that thetotal burden of the anellovirus TTV increases post livertransplantation. Furthermore, an increased prevalence of anelloviridaein pediatric patients with fevers was reported recently.

We next compare the virome composition measured for organ transplantrecipients to the composition observed in healthy individuals, not onimmunosupressants or antivirals (n=9, sequencing data available from aprevious study). Here, we compare the healthy composition to thecomposition measured for organ transplant recipients at the start of thedrug therapy (post-operative day one, n=13), corresponding to a minimaldrug exposure, and to the composition measured for transplant recipientsexposed to high drug levels (well after the transplant procedure,tacrolimus 9 ng/ml. valgancicolvir 600 mg, n=68). We find a similarcomposition of the virome for the healthy reference samples and samplescorresponding to minimal drug exposure (FIG. 2B). However, thecompositions of the healthy reference and minimal drug exposure samplesare distinct from the anelloviridae-dominated composition measured forhigh drug exposure samples.

The tacrolimus-based immunosuppressive therapy is complemented withinduction therapy in the first 3 days post transplant (withanti-thymocyte globulin, daclizumab, or basiliximab) and the patientsfurthermore receive the corticosteroid prednisone throughout the posttransplant therapy. The time-dosage profile for prednisone andtacrolimus are similar: high doses at the onset of the therapy followedby a gradual dose reduction (FIG. 7A-C). The data in FIG. 2A thusreflect the combined effect of prednisone and tacrolimus. An analysis ofthe differential effect of prednisone and valganciclovir on the viromecomposition (FIG. 7E) shows the same trend observed in FIG. 2A: higherprednisone doses lead to a larger representation of anelloviruses.Lastly, we note that a subset of patients was not treated with antiviraldrugs. The data corresponding to this subset of patients allowed us tofurther disentangle the differential effect of the antiviral drugs andthe immunosuppressants on the composition of the virome, as describedbelow.

Partitioning of microbiome diversity. We studied the diversity of thebacterial and viral components of the microbiome. The within-subjectdiversity was lower than the between-subject diversity, both forbacteria and viruses (Bray-Curtis beta diversity, bacterial compositionat phylum level, family and order level viruses, FIG. 2C). Partitioningthe data for patients according to transplant type, heart or lung, orage did not reduce the diversity. Within subjects, the diversity waslower for samples collected within a one month timespan, again both forbacteria and viruses. For viruses but not for bacteria, we find that thediversity is lower when comparing samples collected at a similar drugdosage (tacrolimus level±0.5 ng/ml, valganciclovir±50 mg). Takentogether with the sensitivity of the population averages to drug dosagein FIG. 2A, we thus find that the composition of the virome for patientsthat are subject to the same drug therapy converges to a similar state.

Dynamic response of virome to drug dose changes. A strong temporalresponse of the virome to changes in drug dosage is observed, consistentwith the sensitivity of the virome composition to drug dosage. FIG. 3Ashows the time dependence of the relative genomic abundance of ssDNA anddsDNA viruses (data from all patient groups and samples, n=656). Thefraction of ssDNA viruses expands rapidly during the first months posttransplant followed by the opposite trend after 6 months. FIG. 3B showsthe time-dependent relative composition of the most abundant virusesgrouped at the family and order level and provides more detail on thevirome compositional dynamics (data from all patient groups and samples,n=656). The dsDNA fraction consists of caudovirales, adenoviridae,polyomaviridae and herpesvirales, which together occupy 95% of thevirome in the first week(s) post transplant. ssDNA viruses only occupy5% of the initial virome and mainly consist of members of theanelloviridae family. The fraction occupied by adenoviridae,caudovirales and herpesvirales decreases strongly in the first fewmonths as these virotypes are effectively targeted by the antiviralprophylaxis. In contrast, the relative abundance of anelloviridaeincreases rapidly as these virotypes largely escape targeting by theantiviral drugs and take advantage of the reduced immunocompetence ofthe patients (maximum of 84% during months 4.5-6). Six months after theorgan transplant procedure, the opposite trends are observed, consistentwith the reduction in antiviral and immunosuppressant drugs prescribedby the therapeutic protocol.

Compared to the viral component, the bacterial component of themicrobiome remains relatively stable over time, an observation that ismade at the phylum, order and genus taxonomic levels (FIG. 3C, n=656,and FIG. S3). FIG. 3D shows the within-sample alpha diversity for thebacterial and viral genera as function of time (Shannon entropy, onemonth time periods, 590 bacterial genera, 168 viral genera examined).The diversity of observed viral genera decreases at the onset of thetherapy (1.05±0.5 in month 1 to 0.31±0.33 in months 4-5, p<<10⁻⁶,Mann-Whitney U test), whereas the alpha diversity of bacteria remainsrelatively unchanged during the course of the post transplant therapy(2.2±1.14 in month 1 to 2.6±0.85 in months 4-5, p=0.1, Mann-Whitney Utest).

Increase in total viral load at onset of post transplant therapy. Toobtain insight into the effect of therapeutic drugs on total viral load,we extracted the absolute genomic abundance of all viruses relative tothe number of human genome copies by normalizing the genome coverage ofthe viral targets to the coverage of the human genome. For all patientgroups part of this study an increase in total viral load is observed atthe onset of the therapy (FIG. 4A), regardless of transplant type (heartor lung) or age (adult or pediatric) (change in load, 7.4±3, sigmoidfit, black line). Combined with relative abundance data, the total viralload data reveals a net reduction of the Herpesvirales load and a netincrease in anelloviridae load in the first 3 months post-transplant forpatients that are simultaneously treated with antivirals andimmunosuppressants.

The data thus show a differential effect of the combination ofantivirals and immunosuppressants on different virotypes. The data alsoshow a reduction in total adenoviridae load, indicating thatadenoviridae replication is suppressed by valganciclovir, in agreementwith previous studies. FIG. 4B summarizes data for all transplant types,but the same trends are observed when stratifying according to differenttransplant types: adult heart transplant recipients (n=268, FIG. 9A),adult lung transplant recipients (n=166, FIG. 8B), and the pediatricpatients that are treated with cyclosporine as opposed to tacrolimus(n=99, FIG. 9C).

Not all patients in the study cohort received both antiviral andimmunosuppressant drugs: for transplant cases where both the donor andrecipient do not show evidence of a prior CMV infection in a CMVantibody assay, it is judged that the risks of complications due toantiviral prophylaxis outweigh the potential risk of a newly acquiredCMV infection, and the patients are accordingly not treated withantiviral prophylaxis. These patients are thus solely treated withimmunosuppressants. FIG. 9C shows the time dependent viral load andcomposition of the CMV negative cases (n=75). The net effect ofimmunosuppressant-only therapy is an expansion of all virotypes,including Herpesvirales and adenoviridae. Tapering of immunosuppressionleads to a reduction of the total viral load.

Lower anellovirus burden in patients suffering from a graft rejectionepisode. Given the correlation of the anellovirus burden with the extentof immunosuppression (see FIG. 2A and FIG. 4), and given the associationbetween immunocompetence and the risk of rejection, we asked whether theanellovirus burden can be used for the classification of rejecting andnon-rejecting graft recipients. FIG. 5A shows the anellovirus loadmeasured for rejecting and non-rejecting patients as function of timepost transplant. Here, patients are classified as rejecting in case theysuffer from at least one biopsy-determined moderate or severe rejectionepisode, biopsy grade 2R/3A (in red; 20 patients, 177 data points). Therejection-free patients correspond to patients that are not diagnosed tosuffer from a moderate or severe graft damage throughout their posttransplant course (in blue; biopsy grades<2R/3A, 40 patients, 285 datapoints).

FIG. 5A shows that the anellovirus burden is significantly lower for therejecting individuals at almost every time point. We next directlycompared the anellovirus burden for patients at rejection with theburden measured for patients in the absence of rejection. To account forthe time dependence of the anellovirus load described above (FIG. 5A),we extract the anellovirus load relative to the mean load measured forall samples at the same time point. FIG. 5B shows the time-normalizedload for non-rejecting patients (N=208) compared to the load measuredfor patients suffering from a mild rejection event (biopsy grade 1R,N=102) and patients suffering from a severe rejection episode (biopsygrade 2R/3A, N=22). The figure shows that the time-normalized loads aresignificantly lower for the patients at greater risk of rejection.P-values were calculated by random sampling of the population with agreater amount of measurement points,p=sum(median(Ar_(rej))>median(A_(non-rej))/N, where N=10⁴ and A_(rej)and A_(non-rej) are the relative viral loads for the populations atgreater and lesser risk of rejection and non-rejecting respectively(p=0.011, p=0.0002 and p=0.036).

These observations are in line with a view that the risk of rejectionand the incidence of infection have an opposite association with thepatients' immunocompetence (see inset FIG. 5A). The lower viral loadobserved for rejecting patients is thus indicative of a higher level ofimmunocompetence in this subgroup of patients, even though thesepatients are treated with the same immunosuppressive protocol.Patient-to-patient variability in the sensitivity to suppression ofimmune function is known to occur and the lack of predictability inimmunosuppression is an important risk factor in transplantation. Apresently-used commercial assay for the measurement of immunocompetencehas not been found to be predictive of acute rejection or significantinfections. The development of methods for the direct measurement ofimmunocompetence, that can replace or complement existing assays, willtherefore be important. The total anellovirus load recorded in organtransplant recipients could serve as an alternative marker. FIG. 5Cshows a receiver-operating characteristic and tests the performance ofthe relative anellovirus load in classifying non-rejecting and rejectingpatients (area under the curve=0.72).

We have studied drug-microbiome interactions following solid organtransplantation by sequencing cell-free DNA in the recipients' plasma.The data reveal much about the fundamental structure of the human viromein plasma and how it responds to pharmacological perturbation; they alsoshow the relative insensitivity to immunosuppression of the compositionof the bacterial component of the microbiome. These data are useful inthe design and optimization of post-transplant therapeutic protocols.For example, they show that tapering of antiviral prophylaxis frominitial high doses leads to a resurgence of the herpesvirales fraction.CMV DNA load has previously been shown to predict CMV disease relapseand rejection, raising the question of whether patients would benefitfrom longer-term prophylactic therapy.

The marked expansion in abundance of anelloviridae uponimmunosuppression is also worth further consideration. Anelloviruses areubiquitous in the human population and, although no pathogenicity hasbeen established, anelloviruses are currently under investigation aspotential cofactors in carcinogenesis. The sensitivity of anelloviridaeto immunosuppression makes organ transplantation an ideal setting forthe study of the properties of anelloviridae, particularly in the lightof the increased incidence of cancer seen in transplant recipients. Theobservation of a lower-than-average burden of anelloviruses in patientsthat suffer from a rejection episode is indicative of insufficientimmunosuppression in this subgroup of patients, even though thesepatients were subject to the immunosuppressant levels prescribed perprotocol. This suggests that there would be value in designing assaysthat allow directly gauging the level of a patient's immunocompetence,in addition to measurements of circulating drug levels. The total burdenof anelloviruses identified in a transplant recipient's blood may serveas one such marker of the overall state of immunosuppression of theindividual patient.

High throughput DNA sequencing finds use in the hypothesis-freediagnosis of infections. This approach is of particular relevance in thecontext of transplantation given the fact that infections occurfrequently in transplantation and are difficult to diagnose inimmunocompromised individuals, and given that sequence analysis canadditionally provide information on the graft health through thequantification of donor-derived human DNA circulating in plasma. Inother areas of infectious disease, it may be of value to developsubtractive methods to eliminate the human DNA and enrich for DNA ofviral and microbial origin.

Experimental Procedures

Clinical sample collection: Patients were enrolled at StanfordUniversity Hospital (SUH) or Lucile Packard Children's Hospital (LPCH),and were excluded if they were recipients of multi-organ transplants.This study was approved by the Stanford University Institutional ReviewBoard (protocol # 17666) and enrollment commenced in March 2010. Fordetails on patient recruitment and post transplant treatment of thepatients see the extended experimental procedures section.

Plasma processing and DNA extraction: Plasma was extracted from wholeblood samples within three hours of sample collection, as previouslydescribed (Fan et al., 2008), and stored at −80° C. When required foranalysis, plasma samples were thawed and circulating DNA was immediatelyextracted from 0.5-1 ml plasma using the QIAamp Circulating Nucleic AcidKit (Qiagen).

Sequencing library preparation and sequencing: Sequencing libraries wereprepared from the purified patient plasma DNA using the NEBNext DNALibrary Prep Master Mix Set for Illumina with standard Illumina indexedadapters (purchased from IDT), or using a microfluidics-based automatedlibrary preparation platform (Mondrian ST, Ovation SP Ultralow librarysystem). Libraries were characterized using the Agilent 2100 Bioanalyzer(High sensitivity DNA kit) and quantified by qPCR. Samples were part of26 different sequencing runs and were sequenced over the course of 22months. On average 6 samples were sequenced per lane.

Posttransplant Monitoring and Clinical Sample Collection. This analysisrepresents a substudy of a prospective cohort study funded by theNational Institutes of Health (RC4 AI092673) to study the clinicalutility of a donor-derived cell-free DNA assay for the diagnosis ofacute and chronic rejection and allograft failure after thoracic organtransplantation. Patients were enrolled if they received a heart or lungtransplant at Stanford University Hospital (SUH) or Lucile PackardChildren's Hospital (LPCH), and were excluded if they were recipients ofmultiorgan transplants or if they were followed at centers other thanSUH or LPCH posttransplant. This study was approved by the StanfordUniversity Institutional Review Board (protocol # 17666) and enrollmentcommenced in March 2010.

Details of the Posttransplant Therapeutic Protocol, Adult HeartTransplant Recipients. Posttransplant immunosuppression consisted ofmethylprednisolone 500 mg administered immediately postoperativelyfollowed by 125 mg every 8 hr for three doses. Antithymocyte globulin(rATG) 1 mg/kg was administered on postoperative days 1, 2, and 3.Maintenance immunosuppression consisted of prednisone 20 mg twice dailystarting on postoperative day 1 and tapered to <0.1 mg/kg/day by the6_(th) postoperative month and tapered further if endomyocardialbiopsies showed no evidence of cellular rejection. Tacrolimus wasstarted on postoperative day 1 and dosing was further adjusted tomaintain a level of 10-15 ng/ml during months 0-6, 7-10 ng/ml duringmonths 6-12, and 5-10 ng/ml thereafter. Mycophenolatemofetil was startedat 1,000 mg twice daily on postoperative day 1 and dose adjustments weremade, if required, in response to leukopenia.

All patients received standard CMV (antiviral) prophylaxis consisting ofganciclovir 5 mg/kg IV, adjusted for renal function, every 12 hrstarting on postoperative day 1 unless both donor and recipient were CMVnegative. When able to tolerate oral medications, recipients werestarted on valganciclovir 900 mg twice daily for 2 weeks, then 900 mgdaily until 6 months posttransplant, followed by 450 mg daily until 12months posttransplant, at which point antiviral prophylaxis wasdiscontinued. Valganciclovir dose reductions were made in the setting ofleukopenia. CMV-recipients of a CMV+ allograft also received CMVhyperimmune globulin, 150 mg/kg IV, within 72 hr of transplant, 100mg/kg at posttransplant weeks 2, 4, 6, and 8, and 50 mg/kg at weeks 12and 16 posttransplant.

CMV⁻ recipients of CMV⁻ allografts were not treated with antiviralprophylaxis until May 2012; subsequently, these recipients were treatedwith acyclovir 400 mg twice daily for one year. Antifungal prophylaxisconsisted of itraconazole 300 mg daily for the first 3 monthsposttransplant, and prophylaxis against pneumocystis jiroveci infectionconsisted of trimethoprim/sulfamethoxazole, 80 mg TMP component daily.Prophylaxis against pneumocystis infection was continued indefinitely,and patients intolerant of TMP-SMX were treated with atovaquone,dapsone, or inhaled pentamidine.

All heart transplant recipients were monitored for acute cellularrejection by surveillance endomyocardial biopsies performed at scheduledintervals after transplant: weekly during the first month, biweeklyuntil the 3_(rd) month, monthly until the 6_(th) month, and then atmonths 9, 12, 16, 20, and 24. Biopsies were graded according to theISHLT 2004 revised grading scale (0, 1 R, 2R, 3R) (29). Blood sampleswere collected from heart transplant recipients at the following timepoints posttransplant: weeks 2, 4, and 6; months 2, 2.5, 3, 4, 5, 6, 8,10, 12, 16, 20, and 24. A subset of heart transplant recipients also hadblood samples collected on posttransplant day 1. If blood sampling andendomyocardial biopsies were performed on the same day, care was takento ensure that blood was collected prior to the biopsy procedure.

Pediatric Heart Transplant Recipients. Induction immunosuppressioninitially consisted of daclizumab 1 mg/kg IV every 2 weeks for a totalof 5 doses, and was switched to basiliximab 10-20 mg IV on postoperativedays 0 and 4 beginning in August 2011. Recipients were also treatedimmediately with pulse methylprednisolone 10 mg/kg IV every 8 hr for 3doses, followed by prednisone 0.5 mg/kg twice daily for the first 14days posttransplant; corticosteroids were subsequently tapered offduring the first posttransplant year, in the absence of acute rejection.

Calcineurin inhibition consisted primarily of cyclosporine, with goallevels of 300-350 ng/ml for months 0-3, 275-325 ng/ml for months 4-6,250-300 ng/ml months 7-12, and 200-250 after month 12 posttransplant.Patients intolerant of cyclosporine were treated with tacrolimus.Protocols for prophylaxis against opportunistic infections andsurveillance endomyocardial biopsies were similar to adult hearttransplant recipients.

Lung Transplant Recipients. Posttransplant immunosuppression consistedof methylprednisolone 500-1000 mg administered immediatelypostoperatively followed by 0.5 mg/kg IV twice daily. Basiliximab, 20 mgIV on days 0 and 4, was given for induction immunosuppression.Maintenance immunosuppression consisted of methylprednisolone 0.5 mg/kgIV twice daily on postoperative days 0-3, followed by prednisone 0.5mg/kg daily until day 30, and subsequently tapered every 2-3 months to0.1 mg/kg daily during months 6-12 posttransplant. Tacrolimus wasstarted on postoperative day 0 and dosing was adjusted to maintain alevel of 12-15 ng/ml during months 0-6, 10-15 ng/ml during months 6-12,and 5-10 ng/ml thereafter. Mycophenolatemofetil was initiated at 500 mgtwice daily on postoperative day 0 and dose adjustments were made, ifrequired, in response to leukopenia. Antiviral, antifungal, and PCPprophylaxis were similar to the adult heart transplant cohort.

All lung transplant recipients were monitored for acute cellularrejection by protocol transbronchial biopsies performed at months 1.5,3, 6, 12, 18, and 24 posttransplant. Biopsies were also performed ifindicated for clinical reasons, based on symptoms or pulmonary functiontest results. Blood samples were collected from lung transplantrecipients for study purposes at the following intervals: thrice on day1, twice on day 2, and once on day 3 posttransplant, followed by weeks 1and 2, and months 1.5, 2, 3, 4.5, 6, 9, 12, 18, and 24. Blood sampleswere drawn prior to performance of per-protocol and clinically-indicatedbiopsies.

Workflow for the Identification of Pathogen-Derived Sequences. Exactduplicates were removed using the C-based utility fastq.cpp. Low-qualityreads were removed using the quality filter that is part of the fastxpackage (fastq_quality_filter -Q33 -q21 -p50). The remaining reads weresubsequently aligned using BWA to the human reference genome build hg19(bwaaln -q25). Unmapped reads were collected using samtools (samtoolsview -f4) and low-complexity reads were removed using Seqclean (seqclean-I 40 -c 1). Reads were subsequently aligned to a selection of viral,bacterial and fungal reference genomes and all references in ncbi_fungidownloaded.

FIG. 6A shows the distribution of the genome sizes. The followingparameters were used for the BLAST alignment: reward=1, penalty=_3,word_size=12, gapopen=5, gapextend=2, e-value=10_(_4), perc_identity=90,culling_limit=2. Blast hits with alignment length shorter than 45 wereremoved. For a subset of samples longer reads were available (2 3 100bp, n=55). To test the robustness of the genomic abundance estimates,the length dependence of the composition measurement was examined. Here,reads were trimmed to 40, 50, 65, 80 and 100 bp lengths (fastx_trimmer)and analyzed using the above-described workflow. Here the blast hitswith alignment lengths shorter than 37, 45, 59, 72, and 80 bp wereremoved for the 40, 50, 65, 80 and 100 bp reads respectively. GenomeAbundance Estimation Relative genome abundance estimation was calculatedwith GRAMMy. This tool utilizes the BLAST-derived nucleic acidsequence-similarity data to perform a maximum likelihood estimation ofthe relative abundance of species in the sample. GRAMMy filters hits byBLAST alignment metrics (E-score, alignment length and identity rate)and accounts for the target genome size and the ambiguity of readassignments in assessing the relative abundance of the candidatereference genomes. Grammy was called using following parameters: pythongrammy_rdt.py; python grammy_pre.py -q “40,40,1” input set; pythongrammy_em.py -b 5 -t 0.0001 -n 100 input.mtx; grammy_post.py input.estsetinput.btp.

Custom scripts were used to combine the strain-level abundance estimatesto obtain the abundance at higher taxonomic-level abundances. Here, aminimal taxonomy for the reference database was built using Taxtastic.

Quantification Absolute Viral Load. To quantify the load of infectiousagents in the samples the blast hit results were collected and the besthits selected for each read using a custom script (Bioperl). FIG. 6Bshows the distribution of the number of unique viral, bacterial andfungal blast hits per million unique molecules sequenced. FIG. 6C showsthe number of viral, bacterial and fungal genome copies relative to thenumber of human genome copies present in the sample. The coverage of thegenome of the infectious agent was normalized with respect to the humangenome coverage.

qPCR Validation of Sequencing Results for Selected Viral Targets.Standard qPCR kits for the quantification of Human Herpes Virus 4, 5 and6 and parvovirus (PrimerDesign, genesig) were used to validate thesequencing results for a subset of cell-free DNA samples. qPCR assayswere run on cfDNA extracted from _1 ml of plasma and eluted in a 100 mlTris buffer (50 mM [pH 8.1-8.2]). The plasma extraction and PCRexperiments were performed in different facilities. No-template controlswere ran to verify that the PCR reagents were included in everyexperiment. FIG. 6D compares the relative number of blast hits permillion reads acquired to the concentration of viral genome copies asdetermined using qPCR.

No-Template Control. A no-template control experiment was performed. Asequencing library was generated from nuclease-free water (S01001,Nugen). The library was prepared together with 7 additional samplelibraries (cell-free human DNA) to test for possible sample-to-samplecrosstalk during library preparation. To ensure formation of clusterswith sufficient density on the Illumina flow cell, the sample wassequenced together with a sample unrelated to the study. Whereas thesample unrelated to the study recruited 16 million reads, theno-template control library generated just 15 reads that mapped to twospecies in the reference database, the methanocalcodoccus janaschii (9hits) and Bacillus subtillis (5 hits) genomes. No evidence was found forhuman related sequences, indicating that sample-to-sample contaminationwas low.

Example 2 Clinical Monitoring of the Virome

Using the methods as described in Example, 1 reads that map to the CMVgenome were quantified for each sample. An increased CMV abundance wasobserved in samples that were clinically positive for infection(p=7.10⁻⁹, Mann-Whitney U test, FIG. 10C); the level of CMV-derived DNAin our samples matched clinical reports of CMV with an AUC of 0.91 (FIG.10C). This data indicates that CMV surveillance can be performed inparallel with rejection monitoring using the same sequence data, and ledus to examine whether other viral infections could be similarlymonitored.

We identified well-characterized pathogenic and onco-viruses (FIG. 11A)as well as commensal torque teno viruses (TTVs, alphatorquevirus genus),which is consistent with previous observations of a link betweenimmunosuppression and TTV abundance. The frequency of clinical testingfor these viruses varied considerably, with frequent surveillance of CMV(Human Herpes Virus 5, HHV-5, n=1082 tests in our cohort) relative toother pathogens (FIG. 11A). We evaluated the incidence of infection(number of samples in which a given virus is detected via sequencing)relative to the clinical screening frequency. Although CMV was mostfrequently screened for (335 samples), its incidence as determined bysequencing (detected in 22 samples) was similar to that of otherpathogens that were not routinely screened, including adenovirus andpolyomavirus (clinically tested on four occasions and one occasion,respectively, FIG. 11A).

Adenovirus is a community-acquired respiratory infection that can causegraft loss in lung transplant recipients and poses a particularly highrisk for paediatric patients. Samples were collected from one paediatricpatient (L78, FIG. 11B panel 1) that tested positive for adenovirus.This patient also had the highest adenovirus-derived DNA load across theentire cohort. A sustained adenovirus load was furthermore observed inseveral other adult transplant patients that were not screenedclinically (e.g, L34, FIG. 11B panel 1), as testing is typicallyrestricted to paediatric lung transplant cases.

Polyomavirus is the leading cause of allograft rejection after renaltransplantation but is not routinely included in post-lung transplantsurveillance. We detected polyomavirus in two patients that were nottested for this pathogen (L57 and L15, FIG. 11B panel 2). In both cases,the clinical records indicated persistent renal insufficiency, which mayhave resulted from polyomavirus infection.

In a last example of the benefit of broad and hypothesis-free screeningof infections, we examined a patient that exhibited a high load of humanherpes virus (HHV) 8 (FIG. 11B panel 3), an oncovirus that can causecomplications following solid-organ transplantation. This patient (L58)tested positive for two other herpesviruses (HHV-4a and HHV-5), whichhave the potential to stimulate HHV-8 re-activation. Thoughpost-transplant monitoring for HHV-8 is only recommended in particularclinical circumstances, use of sequencing enables the identification ofthe virus in non-suspect cases that would otherwise go undetected.

Clinical monitoring of the microbiome. In addition to viruses measuredin serum, we also observed correlation between cell-free measurementsand fungal or bacterial infections detected in other body fluids,including Klebsiella Pneumonia infections detected via urine culture(ROC=0.98) and fungal infection detected in BAL. Performance onbacterial and fungal correlations was sensitive both to the infectiontype and body fluid queried. We observed better performance for bodyfluids that have tighter coupling to blood and also observed sensitivityto background signal. For example, the most commonly cultured bacterialinfection (Pseudomonas) was detected in cell-free measurements for over80% of our patient samples, which was in stark contrast to the mostcommonly detected viral pathogenic species (CMV), detected in only 6% ofour patient samples.

This highlights an important distinction between commensal infections(including Pseudomonas), which are part of the normal flora, andnon-commensal infections, which are solely pathogenic and have a lowerbackground signal. This distinction may explain differences insensitivity and specificity measured for commonly cultured commensalinfections (e.g., AUC=0.66 and 0.62 for P. aeruginosa and E. coli,respectively) relative to non-commensals (AUC=0.91 for CMV). In the caseof commensal bacteria, the clinical question is not presence or absencebut rather presence or absence in inappropriate body sites.

In our cohort, we also detected cell-free DNA derived frommicrosporidia, a non-commensal fungus that can cause intestinalinfections in immunosuppressed patients. We measured a sustainedmicrosporidia load in L78 (FIG. 11B panel 4), a patient that exhibitedcanonical symptoms of microsporidiosis. An Adenovirus infection (L78,FIG. 11B panel 1) was the suspected cause, though endoscopy andsigmoidoscopy results were inconclusive and stool samples testednegative for C. diff as well as Adenovirus. Based upon our sequencingdata, the microsporidiosis is the most likely explanation for thepatient's symptoms, as the microsporidia signal measured in this patientis similar to that of 16, a patient from an unrelated cohort that testedpositive for microsporidia (FIG. 11B panel 4).

With more than 10 billion fragments per ml of plasma, circulatingcell-free DNA is an information-rich window into human physiology, withrapidly expanding applications in cancer diagnosis and cancer treatmentmonitoring, genetic prenatal diagnosis, and monitoring of hearttransplant rejection via “genome transplant dynamics” (GTD). In thiswork, we applied the principle of GTD to lung transplantation—aparticularly challenging type of solid organ transplant that is limitedby poor survival rates, as well as an inaccurate and invasive test forallograft rejection.

Because lung transplant recipients with allograft infection and acuterejection may present clinically with similar symptoms, we extended thescope of GTD to infectious disease monitoring. We first demonstrated astrong correlation between clinical test results and cfDNA derived fromCMV—a leading cause of post-transplant graft injury. We further showedthat hypothesis-free infection monitoring revealed numerous un-testedpathogens, including un-diagnosed cases of adenovirus, polyomavirus,HHV-8, and microsporidia in patients who had similar microbial cfDNAlevels compared to patients with positive clinical test results andassociated symptoms. These examples illustrate the benefit of broad,sequencing-based monitoring of infection as opposed to pathogen specifictesting. The approach can be of immediate use as tool that can assist indetermining the occurrence and source of an infection. This may be ofparticular relevance in the context of transplantation, where theincidence of infections is high, where rejection and infection canco-occur, and where the symptoms of infection and rejection aredifficult to discriminate.

While preferred embodiments of the present invention have been shown anddescribed herein, it will be obvious to those skilled in the art thatsuch embodiments are provided by way of example only. Numerousvariations, changes, and substitutions will now occur to those skilledin the art without departing from the invention. It should be understoodthat various alternatives to the embodiments of the invention describedherein may be employed in practicing the invention. It is intended thatthe following claims define the scope of the invention and that methodsand structures within the scope of these claims and their equivalents becovered thereby.

What is claimed is:
 1. A method of detecting a microbe associated with arespiratory infection comprising: (a) providing a sample comprisingnucleic acids, wherein the sample comprising nucleic acids is from asubject with a respiratory infection and is selected from the groupconsisting of blood, plasma, serum, cerebrospinal fluid and synovialfluid; (b) performing high-throughput sequencing of nucleic acids fromthe sample comprising nucleic acids to generate sequence reads thatcomprise sequence reads from at least one microbe associated with therespiratory infection; and (c) detecting the at least one microbeassociated with the respiratory infection.
 2. The method of claim 1,wherein the nucleic acids comprise cell-free DNA.
 3. The method of claim2, further comprising selectively extracting the cell-free DNA from thesample comprising nucleic acids to obtain extracted cell-free DNA,wherein the high-throughput sequencing in (b) is performed on theextracted cell-free DNA.
 4. The method of claim 1, wherein the detectingthe at least one microbe associated with a respiratory infectioncomprises determining a normalized quantity of the least one microbeassociated with a respiratory infection based on the sequence reads fromthe at least one microbe associated with the respiratory infection. 5.The method of claim 1, wherein the nucleic acids comprise cell-free DNA,cell-free RNA, or a combination thereof, and further comprisingattaching adapters to the cell-free DNA, cell-free RNA, or a combinationthereof prior to (b) to produce nucleic acids attached to the adapters,wherein the nucleic acids attached to the adapters comprise nucleicacids from the at least one microbe associated with a respiratoryinfection and wherein (b) comprises sequencing the nucleic acidsattached to the adapters.
 6. The method of claim 1, wherein the nucleicacids comprise DNA.
 7. The method of claim 1, wherein the nucleic acidscomprise RNA.
 8. The method of claim 1, further comprising extractingRNA from the sample to obtain extracted RNA, and wherein thehigh-throughput sequencing in (b) is performed on the extracted RNA. 9.The method of claim 1, further comprising extracting cell-free RNA fromthe sample to obtain extracted cell-free RNA, and wherein thehigh-throughput sequencing in (b) is performed on the extractedcell-free RNA.
 10. The method of claim 1, further comprising identifyingthe at least one microbe at a strain or species level.
 11. The method ofclaim 1, wherein the respiratory infection is a lung infection.
 12. Themethod of claim 1, further comprising, prior to (b), amplifying aplurality of the nucleic acids.
 13. The method of claim 12, wherein theamplifying the plurality of the nucleic acids comprises performingnon-biased amplification on the nucleic acids.
 14. The method of claim12, wherein the amplifying the plurality of the nucleic acids comprisesperforming a PCR reaction with universal primers.
 15. The method ofclaim 12, wherein the amplifying the plurality of the nucleic acidscomprises ligating adapters to the nucleic acids and amplifying withprimers specific for the adapters.
 16. The method of claim 1, whereinthe sample is a plasma sample.
 17. The method of claim 1, wherein the atleast one microbe associated with a respiratory infection comprises avirus.
 18. The method of claim 1, wherein the at least one microbeassociated with a respiratory infection comprises a bacterium.
 19. Themethod of claim 1, wherein the at least one microbe associated with arespiratory infection comprises a fungus.
 20. The method of claim 1,wherein the at least one microbe associated with a respiratory infectioncomprises at least one fungus and at least one bacterium.
 21. The methodof claim 1, wherein the at least one microbe associated with arespiratory infection comprises at least one microbe selected from thegroup consisting of: human respiratory syncytial virus, humanrhinovirus, human influenza virus, human parainfluenza virus,coronavirus, human coronavirus, human SARS coronavirus, humanenterovirus, measles virus, mumps virus, adenovirus, Pneumocystis,Pneumocystis jiroveci, Metapneumovirus, Staphylococcus, Staphylococcusaureus, Multi-drug resistant Staphylococcus aureus (MRSA), candida,Klebsiella pneumonia, Streptococcus pneumoniae, gram positive bacterium,gram negative bacterium, Legionella, Aspergillus, mycobacterium,Mycoplasmatales, Apergillus fumigatus, and Aspergillus niger.
 22. Themethod of claim 1, wherein the sequence reads are short sequence reads.23. The method of claim 1, wherein the subject with the respiratoryinfection is a transplant recipient.
 24. The method of claim 1, furthercomprising quantifying an amount of the at least one microbe associatedwith the respiratory infection based on the sequence reads from the atleast one microbe associated with the respiratory infection anddetecting the respiratory infection in the subject with the respiratoryinfection based on the amount of the at least one microbe associatedwith the respiratory infection.
 25. The method of claim 24, furthercomprising administering an antimicrobial treatment to the subject withthe respiratory infection to treat the respiratory infection, whereinthe respiratory infection is detected based on the amount of the atleast one microbe associated with the respiratory infection.
 26. Themethod of claim 1, further comprising detecting antimicrobial resistancein the at least one microbe associated with the respiratory infectionbased on the sequence reads from the at least one microbe.
 27. Themethod of claim 1, further comprising obtaining a second sample from thesubject with the respiratory infection at a second time point andperforming (b) and (c) on the second sample.
 28. The method of claim 27,wherein the second time point is a time point following administrationof an antimicrobial treatment to the subject with the respiratoryinfection, and wherein the method further comprises detecting a decreasein quantity of the at least one microbe in response to the antimicrobialtreatment.
 29. The method of claim 1, wherein the sequence reads thatcomprise sequence reads from at least one microbe associated with therespiratory infection comprise at least 107 sequence reads.
 30. Themethod of claim 1, wherein the sequence reads that comprise sequencereads from at least one microbe associated with the respiratoryinfection comprise at least 108 sequence reads per sample.